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ABSTRACT 

An analytical technique for the prediction of the effects of 
rigid baffles on the stability of liquid propellant coabustors is 
presented. This analysis employs both two and three dimensional cor~ 
bus'tcr models characterized by concentrated combustion sources at the 
chamber Injector and a constant lloch number nozzle. An elgenfun':tion- 

matching method is used to solve the linearized partial differential 
equations describing the unsteady flow field for both models. Boundary 
layer corrections to this unsteady flow arc used in ■ mechanical energy 
dissipation model to evaluate viscous and turbulence effects within 
the flow. An integrcl stability relationship is then employed to pre- 
dict the decay rate of the oscillations. 

Results of this analysis agree qualitatively with experimental 
observations and show that sufficient dissipation exists to indicate 
that the proper mechanism of baffle damping is a fluid dynamic loss. 

Ihe response of the dissipation model to varying baffle blade length, 
mean flow Mich number, oscillation aiplltade, baffle configuration 


and oscillation mode is examined. 
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Section I 
INTRODUCTION 

A complex sequence of chemical and physical processes takes place 
during combustion and under certain conditions can couple with the 
associated fluid dynamics to produce oscillations in the thermodynamic 
variables and velocity flow field. Random, small amplitude fluctua- 
tions charact irized by the presence of turbulence are typical of this 
flow. A distinct oscillatory behavior, in addition to the turbulence, 
may also be present. These oscillations are organized with a distinct 
frequency and can possess an amplitude which grows with time. This 
type of unsteady behavior is termed combustion instability. 

High frequency instability or resonant combustion produces several 
detrimental effects and is of major concern in liquid propellant rocket 
combustion chamber design. These instabilities have frequencies 
typically between 1000 to 15,000 hertz and have been Treasured at ampli- 
tudes between 10 to 1,000 percent of steady state values. Wall com- 
patability poses a problem with the occurence of this phenomena due to 
increased strecs and heat transfer rates. These effects can then lead 
to rupture and thermal failure of the chamber walls. Other secondary 
problems arise from destruction cf controls and safety devices placed 
internally within the chamber, decreased performance, uncontrolled im- 
pulse and variation in thrust vector. 

As engine development has progressed, combustor designs have re- 
quired Increased performance characteristics and have utilized more 
energetic propellants with Injectors designed to promote more efficient 
combustion. These Influences tend to encourage the occurrence of 
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coijbusMou irstabillty. Instead of altering these cosbustion charac- 
teristics, gachanlcal danping devices have been used to iatprove stabil- 
ity. Two such devices have beer successfully used to suppress 
instability. 

One of these devices is the acoustic liner. It is a series of 
Helmholtz resonators or cir^-saferentiai slots that are placed or ma- 
chined on the periphery of the chamber. Jet losses are responsible 
for the damping that is produced L ' these devices. Experimental veri- 
fication of '.his mechanism has been established^ and a strong theoreti- 
cal basis for design of these devices has been established in several 

2 3 

combustion Instability analysef . ’ However, designers are reluctant 
to use this mechanical damping device because it creates local hot 
spots on the chasd>cr walls and heat transfer becomes an important con- 
sideration. 

The other device that is used is termed a baffle. It is a series 

of blades attached to the Injector surface protr-idlng axially down the 

4 

chamber. This device was first proposed in 1954 and was for some 
time regarded as the panacea re the combustion instability probl^. 
Su'^ficient experimental verification of the stability impr vemer.t of 
combusters with baffles has been produced, however total reliability 
of this device has been limited since a few combustors have failed to 
gain stability improvemeni. with the addition of a baffle. A Lheoretl- 
cai treatment of baffle damping ..s needed to avoid these anonalies and 
ala in denign. Unfortunately, no satisfactory theory exists and design 
of the baffle has remained a black art which utilizes several empirical 
rules that may cr may not be applied effectively in a particular engine 


conf ifvration. 
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Several baffle conf iguratioas have beeu conceived by designers 
(refer to Figure 1) particularly with respect to blade arrahgeaent and 
blade shape. Hcvever, utilizing these designs requires espenslve and 
tim«; ccosualng full scrle tests. Heat transfer aspeccs of these de- 
vic(-z are also of importance but because the baffle is ac internal 
device separable from the ccnabustor vails these ccr.sideraticns aren’t 
cr~<Mcal as far as the structural integrity of the chaaber is concerned. 

Before discussing the theoretical attaspts to model tj't problem, 
the relevant experimcutal observations will be examined. First, ic is 
observed that stability is likely to improve with the addition of the 
baffle. Fur'henaore, an increase in blade length generally increases 
the stability of the chamber. Care with respect to this rule ssust be 
taken since a baffle can be too short or too long." Figure 2 depiccr 
experimental measurements of decay rate for various baffle blade 
lengths.^ This decay rate is representative of the cscillations vi th- 
in a baffled combustor experiencing flow or combustion conditions th.-t 
would produce neutrally stable oscillations in the chamber without 
baffles. Two separate flows, cold flow (flow without combustion or 
mean flow) and hot flow conditions are examined in Figure 2 and show 
similar stability behavior. 

Secondly, it is observed that the addition of a baffle to ? 
chamber depresses the preferred frequencies of the cscillations within 
the chamber. The frequency becomes even smaller with an increase in 
blade Length (refer to Figure 3). In conjunction with these results it 
is noted that baffle contiguration has a minor influence in the problem, 
(several baffle configurations are depicted i.i Figures 2 and 3), pro- 
viding that -he configuration does not coincide with the tangential 
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velocity codecs downstreaa of the baffle. For exaiq>le, an evenly 
spaced thr^e blx-ded baffle has little laflueoce on das|iiag a third 
transverse aode oscillation because wbe velocity nodal lines (in the 
circivaferential direcLinn> are coiixrident vitli the baffle blades and 
the r-^sultins wa - description vithiu the baffle coopartnencs and sain 
chacbe^ are iJ-actical. 

Cold flow acoustic tests conducted at NASA Levis by Vieber^ have 
produced further inforaation about baffle dauplng. Test chafers coo- 
strucced without nonbustlon and nozzle influences indicate the discing 
s'-echaclsz is a fluid dynaalc loss vfaich can be icdep^ident of coabus- 
tlor) and nozzle effects. T.J.S result is also displayed in the sisilar 
suability trends in the cclu flow tests and hot firing tests given in 
Figures 2 and 3. Fron this study it 1st -^Iso concluded that the in- 
crease in surface area and the associate viscous loss (as predicted by 

g 

an idealized theory) is ins*.;ff iclent to account for t- ? energy loss. 

In fact, decay rate measuraaects are an order of magnitude larger than 
predicted by theory. 

A comprehensive theoretical treatment of baffle da&pi^ hac been 

fruitless In previous investigations. Sc*veral analytical EodeVs have 

failed because of oversimplification of the ntot-ies. One suc.l study 

9 

performed by Reardon simplifies the gtjonetry of the pre-bJe® and re- 
quires only one dimensioncl esc illations in the baffle cavities which 
Interact vie:? three aimensionai osclilatlorjs downstrecK of the baffle. 
The influences of coabvstlor., lY^zzle and mean j iow effects ere treated 
and a dsccy rate ralcuiation is made. Frequency prediction with this 
model Is in agieem^nt with experimental observatict, however, results 
from this study inilcate a mechanism for the oaffle damping that is 
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totally 'iependeat upon another loss producing device (l«e. a nozzle) 
in the systea. This conclusion doesn't account for daaping in the 
pure acoustic flow situation which has been observed in Vieber's work 
■entioned earlier. Also the stability predictions are underestiaated 
and do not suov the proper treisis- 

A second aodel that was considered for analytic treataent has been 
suggested by Sirignano and btrahle.^^ This aodel eaploys a distorted 
injector surface (resad>iing a particular baffle coaf iguration) and 
treats a surface perturbation probless. The eq*iation governing the re- 
sulting Clow field are given and an asyrqitotic representation of the 
solution is obtained. However^ since no energy loss or gal- considera- 
tions are aade, stability prediction with this solution is futile. 

The final analysis to be discussed is a aodel oevised by Oberg, 
et. al.^^ This study treats three diaensicnal oscillations everywhere 
in the baffled chaober. An inviscid flow with the influences of cuo- 
centrated combustion, nozzle and mean flow is studied using a varia- 
tional Green's function method. Separate soluti >ns are expanded within 
the baffle compartmenta and matched with the solution downstream of the 
baffle. Pressure and axial velocity are approximately matched at the 
interface connecting these regions and produce the complete solution. 
Owing to the variational method, these matching conditio!:s are approxi- 
mated and in certain cases are grossly represented. However, frequency 
trends are properly predicted ard the acoustic flow fie.d is in agree- 
ment with experimental observation. This solution also predicts a 
strong flow ar i.-d the battle tips which is in agreer.enE with experi- 
mental observation. The major tailure of this theory is its improper 
stability prediction of baffled cliambers. The solution predicts a 
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pressure rise at the Injector end of the chamber and a pressure loss 

12 

at the nozzle. Using the Rayleigh criterion (see also Chapter V) a 
jestabillzlng Influence for the baffles Is then indicated. In other 
words, more energy can be added at the injector and the nozzle ex- 
tracts less energy which results in an unstable flow. 

In revlewlr^ Che foresientioaed theories, it is apparent chat wave 
-ulceration is not the mechanism of baffle damping, falling as indicated 
by the Rayleigh criterion. A second damping mechanism has also been 
considered and a chaise in combustion characteristics. However, 
Wieber's experimental evidence indicates this oannot be Che only mech- 
anism. Concentrating the combustion at Che injector over»phasizes 
the energy input in Che flow. An investigation of distributing the 
combusclon shows that it is insufficient to relax Che combustion energy 
input and baffles still cause an overall destabilizing influence. 

A third possible mechanism is s'uggested from Oberg's work and 
*ieber's experimental res-uits- Since a strong flow near the baffle 
blade tips is indirated in Oberg*s solution, viscous and turbulence 
losses may produce sufficient energy dissipation to overcome the 
driving effect of wave alteration. This fluid dynamic loss Is un- 
accounted for in Wieber’s attenpt to explain his experimental results, 
since it occurs locally at the tips of the baffle blades and departs 
from Idealized theory. In fact, this dissipation mechanism has been 
neglected in all previous theories. 

The incxvislon and estimation of this fluid dynamic loss is the 
purpose of this analysis. An eigenfunction-matching method, which 
parallels Oberg's solution is used to represent the unsteady core flow. 
An exact match, in contrast to Oberg's variational approximation. 
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produces the desired soiutlo:^. Boundary layer corrections In the 
baffle tip region arc added to the solution and applied to a aechanlcal 
energy dlsrlpctlon sodel. 'dssults fros this study show the Inportance 
of this energy loss aochaulsa. 
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Section II 

THBOkY— TWO DIMENSlONAt, CbAMBIR SOUTTION 

Ovlng t<^ the matheaiatical coraplexlty cf this problem, a simplili- 
cation of ccmbustor geometry to a two dimenrlonal rectangular chassber 
is mada. Tils choice of geometry is used purely as a diagnostic tool 
for the Bioie realistic three dimensional flow, and as such predicts 
gross stability trends. However, this model clearly shows the mathe- 
oatical fcrndations of the problaa of interest and application to the 
three diaensioncl chamber is merely a mathematical extension. This 
analysis follows in the next chapter. 

The baffles en^er the problem as discontinuities which rigidly 
protrude axially downstream of the injector end of the chamber. The 
chamber is then split into multiple evenly spaced rectangvlar compart- 
ments which are termlnat'id at tha baffle ends by a single mam chamber 
(refer to Figure 4). 

Combustion and nozzle influences enter the problen. as gain-loss 
boundary conditions. The combustion is assumed to be concentrated at 
the injector face. Support for this assumption is based upon experi- 
mental observation that the majority of the combustion processes are 

13 

completed very near the Injector of the chamber. This model for the 

combustion also overestimates the energy input to the iiow and thus 

14 

represents the worst condition for stability. The unsteady model for 

combustion mass generation used here is assumed to be only pressure de- 

] 5 

pendent according to the Crocco n-T time lag theory. 

On the opposite end of the chamber is a "short”, quasi-steady 
nozzle. Due to the restrictive nature of the flow within this nozzle 
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a constant Mach number condition exists at its entrance. The choice 
of this loss boundary condition has a secondary Importance in this 
problem because the nozzle produces a minor damping influe:ice. 

Periodic oscillations are treated for a thetikslly and calorically 
perfect gaseous flow. The cjncentrated combustion assumption permits 
the gas dynamic flow field to be represented as a single constituent, 
product gas with no heat transfer or diffusion processes taking place. 
The core flow within tT<e chamber is characterized by a constant Mach 
number steady flow and is devoid of molecular viscosity and turbulence 
effects. Consequently it is consistent to assume a vexocity potential 
exists for the core main flow. Corrections to these assumptions are 
made by making boundary layer aujustments at the appropriate surfaces. 

One final assumption is nade with regard to this sol«itiou. 

Entropy variations a-e neglected in this analysis. This assumption is 
consistent with the small overall influence they produce on the 
problem. 

Before mathematically describing the preceeding flow, a non- 
dimensionallzation of the thermodynamic variables and the velocity 
field with respect to the steady state values is made. Because of the 
concentrated combustion assumption, the steady state thermodynamic 
variables and gas velocity are spatia’ly independent and are represent- 
ed as constants. 

The nondimenslcnal conservation equa'.lons governing the flow are 
given in the following relationships. The conservatloii of mass has the 
form: 

^ + o^*q * 0 

n#' • 
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where p and q are the densltv and gas velocity. These quantities are 
nondlaenslonallzed as: 


P * P /p 

-*■ .»•* — -i 

q = q /a 


-* 


where a is the dlEenslonal steady state speed of sound. The Independ- 


ent variables z, y and t are nondlmeasionallzed as: 

*, * 

y = y /R 

*. * 

z “ z /R 


t = 

* 

In this two dimensional chamber R 
* 

along y . 

The momentum equation is: 


it—* 1C 


t a /R 

refers to a characteristic length 


P 




0 


where P is the aondlmcnsional pressure given as P = F /P and y is the 
ratio of the specific heats. 

The final relationship is the homeatropic condition: 



The velocity potential assumntion allows the velocity field to 
have the following representation: 

- 5 ♦ 

The state variables a e then represented as power series expansions 
in an amplitude parameter (e), i.e. 

9 = 5 + e <4>' + O(e^) 

P = 1 + eP' + O(e^) 
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With these expansions a first order linearization of the conser- 
vation equations is made which yields the following equations: 


P' - 


p' « yp' 


dt* 3z3t 3z^ 


- Tfl*' 


I 3t 


J 


( 1 ) 

( 2 ) 


Since standing wave solutions are examined in this analysis it is 
consistent to assume exponential time dependence of Che per cur bat ions. 
Therefore the perturbed pressure and velocity potential have the forms: 

P' = PWe^'*^ and 4»' = (^(R)e^“*^ 

where u - + iX is the complex frequency aiid X is the decay rate. 

Substitution of this time dependence transfoms equations 1 and 2 into 
the following forms: 

V^4) + = 2iw M 1^ + (3) 

aZ az 

P = - Y j i(4 + M II] (4) 

The gain-loss boundary conditions at the injector and nozzle 
entrance surfaces are fomwlated, respectively, with the aid of Crocco’s 
n-T time lag theory and the 'short" nozzle approximation. The combus- 
tion boundary condition is mathematically expressed as: 

m = m n -j ? (t) - P (t - t) 
or 


Ml 

3z j 


z»0 



,, -IWT. 

M n(l-e ) 


z=0 


(5) 
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where n Is the Interactioa index « (a measure of the amplitude depend- 
ence of the mass generation on the pressu'e), and T is the sensitive 

time lag, (a time phasing of the mass generation with the pressure). 

— 18 
Typical values for n and T take the following ranges: 

-* 

0.14 ^ X < 0.20 mllllsecs 
0.6 ^ n £ 0.8 

The "short" nozzle (constant entrance Mach number) approximation 
allows a loss boundary condition which is expressed mathematically as: 


“ p| (6) 

z=L 'z*L 

On the remaining surfaces of the chamber and on the baffle blade 
surfaces hard wall boundary conditions are used. This is expressed by 
a zero normal component of velocity. 

With the partial differential equation for the first order veloc- 
ity potential and the linearized boundai'^ conditions a solution is tJien 
obtained. Because or the dlsconti"uous geometry of the problan, a 
separation of variables solution can not be directly obtained and a 
more sophisticated method is necessary. This method calls for separate 
solutions in the baffle c ivlties and the main chamber. A matching of 
these solutions is made at an artifical Interface between these regions 
thus producing the complete solution. 

The solution within the baffle cavities (0 < z < z ) is found by 

- — B 

separation of variables and utilizes the injector and hard wall 
boundary conditions. This solutio.i takes the form; 


3i| 

3z| 


'P’ = 


w 

E 

m=0 


A^ COSmJlNy 


D 


e^^^’B^B + Cg e^®’’B^Bj 


(7) 
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lAere N is the total number of baffle cavities and^U is an integer 
ifliicb defines a particular baffle coiqsartnuit for £ 7 5. W/K 

(refer to Figure 4). The constants * ^ 2,3 ^3 defined 

by the differential equation and the boundary conditions. (These con- 
stants are explicitly defined in the Appendix.) 

Hithln the main chamber (z < z < L) the perturbed velocity 

B — ~ 

potential tfkes the form 




Oft 

^0 


COSnory 




( 8 . 


where B , B and C a_c constants defined by the differential 
i,c 2,C c 

equation -.nd toe nozzle boundary condition. (These constants are given 
In the Appendix.) 

The complete solution is then obtained by determining the proper 
set of I and 

potential at the main chamber-baffle compartment interface produce the 
conditions necessary for the specification of the eigenfunction 
coefficients. 

With the aid of the orthogonality properties of the series the 
following matching equations are obtained: 


{“»}• The continuity of axial velocity and velocity 
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and 



14 


fl_»> 


iB + 13, C 
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^lB.^(z-L)^ (z-M r 
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/“■ ■a <^B \ 
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+ Cjj e^®2,B*B 

«/» 

M / COSn'Try OOSmirNy dy 

— - 

, J COS^'ny dy 


(10) 


r.':* ,tion 9 is the representation of the matching of velocity po- 
tential and Equation 10 is the matching of axial velocity at z=z^. It 
is recogn:>ed that the solution to this problem satisfies a homogeneous 
differential equation with homogeneous boundary conditions and as such 
poses an eigenvalue problem. Since the amplitude is arbitrary in this 
solution, a normalization to a particular mode within the main chamber 
is made. This gives an additional relationship that is used to com- 
pute the eigenvalue (frequency). Mathematically this is expressed as: 


B.= 

ffl 
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or 




,c^^B + IR. C e^®2 c 
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GOSntry CO&nrNy dy 
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( 11 ) 


OOS*Dny dy 


where a refers no the domlcatiug transverse >node in the main chamber. 

A successive approxlmatloc technique is used to solve Equationa 9> 
10 and 11 . The first approximation chosen for this method uses the un- 
baffled chamber velocity potential soliition, i.e. ® ^ * h’ith 

this approximation a calculation of the baffle compartment coefficients 
^A^l-is made using Equation 9. These coef ricients, in turn, are used 
to recalculate the main chamber coefficients fnom Equation 10 and 

the eigenvalue Equation 11. The procedure is Chen repeated until con- 
vergence is obtained. This iteration scheme converges very qt ickly 
and produces frequency predictions which have less than 5% differences 
after approximately 5 Iterations. 

In investigating the convergence of this solution, thr tching 
relationships are becked by examining the velocity potential and axial 

velocity predictions at the interface (z=z ), Figures .nd 6 show 

B 

these plots for a two and three compartment baffle ronflguratio . With 
the use of Cesaro summation of the series expansions, reasonable agree- 
ment is obtained to assure correctness of the mathematics. These 

figures also indicate large velc ties at the baffle blade tips (z=z ). 

B 

Ac ihese regions the eigenfunction expansions fall to accurately 
represent the flow field. 


It is then necessary to characterize the velocity field near Che 
baffle blade t^ps. To treat this problem a polar coordinate system is 
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set at the blade tips and an expansion of the velocity potential is 
c'btaieed (refer tc Figure 7 for ebe coordinate systea). 

The partial dift'erectlal equation describing the unsteady flow 
Che tips , £q«^tlon X is cransfotMd into the following equation; 






M*# = 2iii*l 


r Sinx 3d 
I C 3o 




The blade boundary conditions are represented £s 


^“U-0,2ir 


0 


If a proper ordering of the solution is aade with respect to an 
asyeptotic fora, l.e. 

^ = ?^A(a) + k 

where 0 ^ ^ « 1 and s ^ 0, it is found that the mean flow corrections 
(right hand side of above partial differential equation) are of the 
orders 0(CM) and 0(M^). Also the <i>^$ tern on the opposite side of the 
equation is a term of 0(<;^). These terms are very .’imall and are 
neglected to produce the asymptotic solutions: 

(j> “ a C** COSa/2 + k + 0(C) (12) 

C « 1 

This solution indicates tl sing- -r behavior of the velocity 
field, i.e. 3<f/3c and 1/C 3ct-/3a *; C ■*<> • The constants a and k are 
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detenlaed bv aacching this asyaptotic solution at soae C vlth the 

c 

outer eigenfunction expansions. Proper choice of ^ is aade so that 

c 

the remainder teras of Equation 12 are negligible and a region of 
proper oatching of the outer expansions is realized. 

Since only two aatch points are needed for this inner solution, 
the asyantotlc expansion and the outer series expansions for the range 
eacircling the blade tip are checked for consistency. Figure 8 shows 
this co^iarlsbn. The velocity potential is well represented by this 
as 3 TipC>tic expansion particularly for v/2 £ a ^ 3»/2 . The differences 
at the otJter yarts of the range are iaaediately accounted for by recog- 
I ling that Che series exp^^sions have outer boundary conditions. 

This, however, is not of •sajoi concern (as will be explained later) 
since the region of Intert^t Is it/2 ^ tx 3x/2' - 

with this representation of the strong flow near the baffle tips, 
a calculation of energy dissipation due to molecular viscosity and 
turbulence is then made. Boundary layer corrections on the velocity 
are necessary. This correction for a laminar, pezicdic flow is given 

U(^,t)=U(rpS,t)(l-e''^'’) fl3) 

where UCn*5,t) is the periodic outer flow transverse velocity and n is 

the normal component to che boundary surface (refer to Figure 9). 

An e.'tinat€ of the meciianlcal energy dissipation within the bound- 

20 

ary layer voJuae Is then calculate! fiom the following integral: 


dis 



dndb 


,//¥ 


U*(n=6,t) ds 


(U) 
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Steady flov correction are neglected in this calculation since 
these corrections are an order aignltode saaller than the unsteady 
flow velocities. This is also consistent with the neglect of dissipa- 
tion on the chasber walls. 

A tiiae average of the above relacioashlp will then physically 
represent an average quantity of aechaaical energy that is transfomed 
irresversible into heat. 

With the presence of conbustion a highly turbulent flow situation 

atist be realized and the existence of turbulence prodt^es aore energy 
21 

dissipation. To account fox. this dissipation, the Boussinesq approx- 
Inatlcc is retained which uses a stress-strain law fo- the tiae aver- 
aged turbulent flov. A "turbulent viscosity" which is function of 

the local flew conditions is then necessary for the aodel. Many rela- 

7.2 

tienships exist for this partaieter, each having liihxted applicability. 

These relationships require a steady flow (apart from the turbulent 

fluct’iatlons) far fron the boundary surface. Models incorporating 

unsteady outer flow are nonexistent and Bust be created froa existing 

steady flow theories. An effective viscosity aodel created by 
23 

Spalding Ic used for this analysis. This aodel is chosen because of 

its simplicity and Its qualitative accuracy with respect to other coa- 

24 

bustion flov probleris- ' This model is represented as: 


* 

^eff p 


* 2/3 


• ^ 


*2 


• *2 

o U 1 
0 0 ^ 


wher*. the F and 0 subscript refer to the fuel and oxidizer quantities. 


In order to be consistent with the single component gas assumption of 


this problem, the fuel and oxldizec veloc: ties are assumed to be uni- 

* * 

form with the product gas velocity, i.e. = U then the 
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equatioQ reduces to: 


eff 




■ *1 

'where |U ] is a r.a.s. value of speed In the entire turbulent field. 

An Incorporation ^f the periodic flow velocity ln;:o the aodel i:i 
then aade by tlae averaging the r.«.s. jO j to give the final form of 


Geoiaetrlcal corrections to this e^juatlon are inflected since they 
have ft weak dependence in the aodel. Spalding suggests a proportion- 
ality constant (C ,) of rhe ordei. 0(0.05). 
t.uri> 

This Bodel is then used to calc'ilate the turbulent dissipation. 

By assuning that the turbulnct velocity profiles are sinllar to the 

laainar predictions given by Equation 13 » th? following integral r 2 la- 

25 

tlon for the dissipation Is obtained: 



(16) 


The isportance of the exactness of the turbulent velocity profile is 
secondary since a global, integral quantity is evaluated. However, it 
is experiaentally observed that turbulent profiles are steeper in shape 
than the laminar flow profiles and consequently this dissipation cal- 
culation could be underestimated. 

Cne final correction is necessarj' for the dissipation calculation. 
A physically impossible infinitely thin baffle blade will create an in- 
finite amount of energy loss because o.'^ the singular behavior of the 
volccity at the tip, A baffle blade of finite thickness will therefore 
be used in this problem. 
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To correct for this thicknese a neighboring strealine is nsed to 
represent the baffle surfaces. This eliwiaates a reworking of the 
solution to correct for the baffle shape since the normal component of 
velocity vanishes along a streamline. Figure 10 shows a streamline 
plot of the fl<w within a baffled chamber with no mean flow. Near the 
tips of the blades the streaallnes are veJJ. riqiresented by 
^ * C S^aqt/2 particularly for ? « 1 and 's/2 ^ 5. 3v/2 . This 

defines tbr geometry of the blade tips. This streamfunctlon 
re^^iued to describe rest of the baffle blade surface but because 
the velocity decreases substantially away from the blade tips this sur- 
face description la of secondary importance in the dissipation calcula- 
tion. Ha thematically this surface is represented as: 

^ Slno/i « (T/4)^ 
where T is the blade thickness. 

A calculation of the dp loss is now available using Equations 12, 
15 and 16. Rather '.'nan correcting boundary conditions to account for 
this dissipation, as is done in aco;istic theory, a nore direct method 
of stability prediction is applied. An integral time average of the 
energy equation, derived by Cantrell and Hart,^^ Is used to estimate 
the global stability of the flow within the combustor. Stability 
behavior (a calculation of decay rate of the perturbations) is examined 
by accounting for the energy inputs or extractions at the various sur- 
faces of the chamber. Mathematically this relationsjjip, correct to 
is stated as: 
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+ 2 




(17) 


wht>re X is the decay rate that is of the order 0(e). 

Evaluating the r.h.s. over the baffle surfaces results in a tera 
representing the aechanlcal energy extracted at the surface. This can 
directly be equated to the dissipation integral Equation 16. By 
applying the appropriate boundary conditions the Cantrell and Hart 
Integral relationship has the final fora: 
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(18) 


dV 


I 

I 


The first tens in the moBeracor represents the energy added to the 
unsteady flow by the combustion, the second term the energy cjctracted 
by the nozzle, and the third term the energy loss created by the strong 
flow surrounding the baffle blade tips. 

Stability calculations using this mathematical analysis are coded 
in Fortran and evaluated using a CDC 6400 computer. Inputs to the 
program Include combustor geometry, baffle configuration, mea flow 
Mach number and oscillation mode character. Some aspects of ti,a com- 
putation problem are discussed later. 
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Section III 


THEORT—THREE DDfENSIOHAL CHAMBER 


A nore realistic description of the prob^eo of Interest is made 
by discarding the two dimensional geraetry assumption and treating a 
three dimensional cylindrical combustion chamber. Tba baffles enter 
the geometry in a similar manner by dividing the injector end of the 
chamber into equal angle sector compartments (refer to Figure 11) . 

Since the same partis differential equation and boundary con- 
ditions are applied to this problem, the extentlon to the three 
dimensional problem reqi.lres doubly infinite eigenfunction expansions 
in the main c!iamber and baffle compairtraents. These expansions are as 
follows: 


m=L jc=l 




(19) 


and 


00 cc 


■ T.T. 

m=0 1=1 
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( 20 ) 


where 

«>SmNe 

2 2 


2]lC; -1) 1 0 ^ 2tt)' 
N n” 


Within the main chamber two types of solution are possible. The 
standing wave solution is made by specif ing: 


(r.e) = C0Sm9 J (X^ r) 
£m m m 



23 


The traveliag wave solution is given by: 




The constants B, „ B, _ C_ B, B and C are given explicitly 
i,B, 2,B, B, »,c, 2.C c " 

in the Appendix. 

The matching relationships and eigenvalue equation that are 

necessary in obtaining a complete solution take the following forms: 
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The eigenvalue equation is represented as: 


B; = 1 
tm 


or 
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where o and £ respectively specifies the dominating transverse and 
radial modes in the main chaiid>er. 

An initial approximation to the preceeding equation sets is sim- 
ilar to the two dimensional chamber problem and is given by 

®im ' 

The same iteration scheme described in the previous chapter is then 
used to obtain higher approximations. 

Figure 12 shows a comparison of the matching of the velocity po- 
tential and the axial velocity at the baffle ccmpartment-maln chamber 
Interface. Reasonable agreement of the matchings js predicted by the 
baffle compartment solution as compared to the main chamber solution 
indicates the correctness of the solution method. Also Indicated is 
the singular behavior of the velocity field at > baffle tips. 

An asymptotic expansion is used to describe the velocity potential 
near the tips with the allowance of radial dependence of the parameters 
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a and k lu Equation 12. With this solution the energy dissipation 
given by Equation 16 is used with the Cantrell-Hart stability Equation 
18 to predict decay rate. 

A Fortran program is used to analyze the results of the calcula- 
tions with the inputs of combustor geometry, baffle configuration, mean 
flow Mach number, and main chaisber mcde character. Cooputatlou is more 
lengthy in the three dimensional problem and requires 20 times the 
computer time required for the two dimensional chamber problem. The 
next section gives some details concerning the computational problems 
Involved . 
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Section IV 
C(»fPUTER S(n.UTIOIi 

Several aspects of the craputer evaluation of the preceeding 
solutions aerit further explanation. First, convergence of the solu- 
tion with respect to approxinating the infinite series expansions and 
the iteration method used to calculate these expansions is discussed. 

A major concern in this solution is the accuracy of the velocity po- 
tential representation. 

In the two dimensional chanber solution, the number of tema used 
to approximate the infinite series is not a critical problem because 
the solution is represented as a single series expansion that requires 
relatively little computing time. Topically a 30 term expansion in the 
baffle compartments and main chamber produced sufficient agre&aent of 
the match of the velocity potential at the baffle compartment-main 
chamber Interface. This U'jmber of terms is also sufficient for the 
match of the baffle compartment series solution to the asymptotic solu- 
tion near the blade tips. Increasing the number of terms has only a 
slight effect on the overall model predictions, decreasing the number 
of terms causes only a modest reduction in computer time required. 

Compared with the two dimensional chamber solution, the three 
dimensional cylindrical chcmber solution is more complex. The compli- 
cations of calculating the solution arise from the double infinite 
Fourler-Bessel expansions by which the solution is represented. The 
number of terms used to approximate these infinite expansions Is very 
critical and drastically Influences computing time and storage. 

Improved convergence and the elimination of Gibbs phenomena has been 
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attained using Cesaro sun&ation. The use of this technique allows a 
reduction in the niuaber of t<irms retailed in the expansions and com- 
puting time is consequently reduced. Figure 13 shows a compirison of 
the baffle compartment and main cbanber series solutions of the axial 
velocity at the baffle-main chamber interface (z*z^). These 1C term 
expansions (susaaed without the Cesaro technique) are grossly matched. 
Figure 14 shows the improvtsaent of the matching with cbe appi icatiou of 
Ce^rc summation to the series solutions. 

The second aspect cf convergence Is concerned with the iteration 
schoe used to solve the matching relationships and the eigenvalue 
equation. This is not a major problem because the sucreseive 
approximation technique converges very quickly. Typically after 5 
iterations the frequency predictions differ between successive Itera- 
tions by less than 5%. However, rhe eigenvalue equation h's multiple 
solutions so care must be taken In choosing the proper initial 
frequency. 

lir arriving at ths solutions several Integral quantities . "e 
numerically evaluated using an improved quadrature method. Eli’ven 
quadrature modes evaluate the dissipation integral with sufficient 
accuracy. \nother Integrant antity is eiraiuated using eleven cuadra- 
ture nodes and appears in the cylindrical chamber matching relatiop- 
shi"s : 

f ■'os 

Jo " ' 

A matrix of these Integrals is evaluated numerically and represents a 
major part of the computing and storage of the cor-puter program. For 
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large values of j^/2 d' quadrature inodes must be used 

to assure proper convergence of the solution. 

Typical computation times for a particular baffle configuration 
(and consequently one and X) are 60 seconds for the two dimensional 
chamber solution with 30 term expansions and 1000 seconds for the three 
dimensional chamber solution which uses iCxlO Fourler-Bessel expansion 
for the solution. Better optlmlk:atlon and storage could possibly 
reduce these computing requirements. 
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Section V 
RESULTS 

Previous investigations base^< on wave alteration as the possible 
mechanism for the damping produced by baffles have failed to correcDy 
predict stability trends. These studies have all neglected viscous 
and turbulence effects In the flow. The incorporation of these effects 
is of major importance in this proble i and is necessary to describe the 
stability behavior of baffled combustors. 

Before discussing results fre the present model which uses me- 
chanical energy dissipation as a damping mechanism, the shortcomlrtgs of 
previous analytical attempts will be clarified, A weli established 
criterion stated by Lord Rayleigh clearly indicates why wave a.teration 
may cause a destabilizing behavior In baffled combuscorf The sense of 
the criterion is giver ty Rayleigh's statement, ”lf heat be given to 

the air at the moment of greatest condensation, or be taken from it at 

12 

the moment of greatest rarefraction the vibration is encouraged." A 

mathematical formulation cf this staremeut has been made by Cantrell 

and Hart and is given by Equation 57 in Section II. Applying the 

appropriate boundary conditions on the injector surface, S and on 

the nozzle entrance surface, S , a! lows an integral inaqualitv which 

noz 

can be used to predict the global stability of the perturbed contined 

flow. By referencing the combustion parameter n and T to those of an 

unbaffled neutrally stable combustor this inequality is reduced to a 

function of only the perturbed pressure distributions on the surfaces 

S, . and S . (In this reference chamber oscillations neither grow 
ii.j noz 

•• dei ay with time because energy is aeded or extracted at the same 
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pressure level.) Formalin, the inequality requires that for unstable 


oscillations: 





In other vords, if aora energy is added due to the conbustlon than is 
extracted by the nozzle the flow is unstable. The equality sign lndi~ 
cates neutral stability behavior of the flea. 

Figure 15 shows the pressure distribution for a baffled and un- 
baffled coad>usLor at the nozzle (z^L) and injector (z=0) ends of the 
chanber. Without a baffle the conbustor is neutrally stable and shows 
dentical pressure distributions over injector and nozzle surfaces. 
Plac»ent of a baffle in a chamber results in a pressure amplitude rise 
at the injector surface and a pressure amplitude loss at the nozzle. 
Consequently, the flow becomes unstable. 

It is this pressure wave alteration which has resulted in earlier 
investigatlcns predicting a destabilizing influence for baffles. Or.Iv 
with Che consideration of mechanical energy dissipation will the sta- 
bilizing Influence of baffles become apparent. Seme results for the 
present model which includes this dissipative influence will now be 
discussed. 


Two Dimensional Cha mber Results 

Linear stability predictions are examined first in a two dimen- 
sional combustor modeled as a rectangular chamber with a leagth L>°1.5. 
First transverse mode oscillations are stvtdied and are assumed to dom- 
inate the solution within th_ main chamber. (First transverse mode 
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refers to a standing oscillation in the j direction that has a half- 
cycle between y=0 and y*l.) The paraaeters n and x required for the 
unsteady co^ustion aodel are assigned the values c»(y^1)/4y and T=l. 
These choices are made so that a neutrally stable unbaffled coisbusvor 
is referenced (i.e. with no baffles the chaaher is aeutrally stable). 

For the dissipation aodal a value for the turbulent coefficient 

e 

23 

C^urb ~ 0-034 is used and is tne value suggested by Spalding. .\lso 
a nondiaenslonal thickness which is typical of actual baffle configura- 
rlons of T*0.05 is chosen for the baffle blades. 

Before examining the stability trends, the effects of blade length 
on the normalized frequency, (frequency of the baffled chacber/f requoi- 
cy of the unbaffled chaid>er), are examined. Figure 16 shews that the 
frequency decreases when a baffle is added to a chanber. Furtheraore, 
an Increase in blade length further depresses the frequency. This 
urediction is in agreesaent with the experioental data provided by 
Aerojet-General. Mean flow corrections to this frequency as predicted 
by the laodel are shown in Figure 17 It is seen that the normalized 
frequency is reduced with an increase in mean flow Mach number. 

The principal result of these calculations is the prediction of 
combustor stability. Decay rates are calculated for various baffle 
blade lengths and chamber conditions. Decay in decloel/cycle is de- 
fined as follows: 

Decay In declbels/cyclc = 20 log^^J^ - p -r 
Wltl. che exponential time dependence of the oscillations In this prob- 
lem this definition reduces to: 

Decay in decibols/cycle = 54.575 
where X is the decay rate and is the frequency. 
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Figure 18 shows decay predictions for a two c<mpirtment baffle 
configuration in a two dlnensfonal cowbustor with nc wean flow (pure 
acoustic flow). Two wave aaplltudes (c) are exaKlned and are shown to 
bracket the qualitative stability trends as reported lu experlsHental 
observations (again froai Aerojet-General data). With an increase in 
anplltude the coabustor stability is iaproved. This isproveaent is 
attributed to toe increased kinetic energy of the flow which increases 
the turbulent viscosity near the blade tips. 

With the addition of coubustion and nozzle Influences la the flow, 
the flow field aust contain aean flow corrections. These influences 
increase the energy densities of the flow and consequently increase the 
sechanical energy dissipation. Bowever, the addition of co«d>ustloi 
creases the energy input intc. the flow and lessens t!ie global sLability 
of the flow. Figure 19 shows the stability trends for a c<^>ustor ex- 
periencing oscillations with waplltude C =* 0.1 for various aean flow 
Mach miners. Even trlth the conbusticn input, sufriicleot dissipation 
exists to stabilize the chamber providing that the baffle Mades are 
long enough. The effectiveness of increasing the blade length of the 
baffle Is reduces as the aean flow Mach number is increased, providing 
the wave ampl:*tude stays constant. Beyoca a blade length Zg > 0.1 the 
nozzle virtually loses all its damping ability because the pressure 
oscillations are distorted so that a pressure node at the nozzle 
entrance occurs. Consequently, only the baffle tip dlsstpaclon reo^alns 
to counteract the energy input at che injector. It can Ve seen in 
Figure 19 ti»at ^or large Mach nuroters a brffle can destabilize a coc- 
bustor because the increased energy due to the combustion overpowers the 
energy loss at the blade tips. This may be an explinatlon for che 
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exper taental cbserv&tion that baffles can be too short and actually 
destabilize the coobustor. However, it aus*' be recognized that the 
concentrated coirimstion aodel overpredicts the energy input to the cos^ 
bustor and in doing so may be responsible for the destabilizing be- 
havior . Mso T/ith an increase in coadnistion input into the chaaber an 
increase in wave aaplitude (e) can be expected. This Mach nuaber de- 
pendence in c is neglected in the results of Figure 19 and consequently 
the dissipation is underestimated for large Mach nuad>ers. 

Stability trends for several oaffle configurations are shown in 
Figure 20. Baffle configuration is seen to have a secondary influence 
in the problem. Adding aore coapartaents to a configuration decreases 
the tip velocities and also increases the surface area. These in- 
fluences counterbalance each other and produce small changes in sta- 
bility character of baffle configuration. There is an advantage of 
large compartaant configurr.tiors, however, since jore modes can be 
effectively damped vlth aacy coapartnents. For example, a two compart- 
ment baffl^i is only effective in damping odd numbered transverse modes 
whereas a five con*partr«ent baffle is effective in damping all modes ex- 
cept those which are multiples of five. 

Figure 21 shows the effectiveness of a two corapartaent baffle con- 
tiguration in damping the first anc chlra transverse mode in tne main 
cVamber (wave amplitude e = 0.1). The third transverse mode character- 
izes a standing oscillation which has three half-cycles between y =• C 
and y = 1. Figure 21 displays the fact that a baffle which damps the 
tirst transverse mode also damps the third transverse mode. The 
decibel rating of decay/cycle is misleacina in this figure because the 
frequency scales are different. (The frequency of the third transverse 
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oode is about three tines that of the first transverse ande.) A aore 
representative stability plot Is giv^ in ?i<rsxe 22 and is a plot of 
decay race (X) for various blvle lengths. This figure indicates the 
sinilar behavior of the two cn;v>a. 

Three Dlnensional C^»a?^er Results 

Of greater- practical interest are the stability predictions for 

the three diae«isional cylindrical chaafcer' 9k>del used in this work. A 

parciculrr cuaid>er with a length to ra/^ius ratio o? L/R * i.5 is 

exa*^.ned. three bladed, evenly spaced baffle configcratlon with 

blade thickness T ~ 0.03 is Input into the dlsslpatlcn nodel for the 

rcrbulent viscosity, ihree diXiensional first transverse nude oscUla-' 

tions in the aaln chanber are studied with an ucsteady co^Mistlon input 

which use n “ (v1-1)/4y and T = ir/X;- . Again a neutrally stable un- 

xa 

baffled coabustor is referencc^d with these paraneters. 

The frequency predictions of the three diaensional solution par- 
allel those of the two dineosioacl solution aixl are in agreoent with 
exporinental data. These results are shown in Figure 23- 

Decay rate predictions are shown i-’. Figure 24 for a coabustor with 
no aean flow aed no cciuburtlon or nozzle Influences- Two wave anpll- 
fjdes, (e = 0.1 and e = 0.2} are snown in this tlgure and Indicate the 
stabilizing behavior of baffles, because of the additional redial de- 
pendence in the three dimensional chember solution the baffle tip 
velocities are not as large as those p-^edicted .^n the two dimensional 
chamber solution. Consequently, the dissipation is overestimated in the 
two dimensional chamber solution relative to the more realistic three 
dimensional chamber case. 
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Figure 25 shows the stability prediction, of a cylirdrical cc<n- 
bustor which has a mean flow (M = 0.1), ci ;&bustion and nozzle influ- 
ences and is experiencing oscillations with amplitude e = 0.1. Two 
types of oscillation are possible for this combustor: standing and 

traveling wave oscillations. Without dissipation Influences, it is 
seru that these solutions produce the same destabilizing Influence for 
Uie baffle (as predicted by the Rayleigh criterion). With the Inclu- 
sion of wchanlcal energy dissipation in the model, these results are 
reversed and show a stabilizing Influence for the baffle. It is seen 
■that the traveling wave solution Is c»st affected by the presence of 
the baffle and produces decay rates that are four times those of the 
standing wave solutio''. This is a critical result because the travel- 
ing wave is most common and is the most destructive. It is also 
apparent tbat the phasing between the oscillations In the main chamber 
and the standing wave oscillations in the baffle cavities produce dif- 
ferent stability results. This observation is also made in Vieber's 
experimental results.^ 

Mean flow corrections to decay rate are shown In Figure 26. These 
predictions agree with the two dimension-*! chamber solutions and Indi- 
cate that a particular baffle blade length becomes less effective with 
an increase In mean flow Mach number, providing c Is constant. Again 
the Mach number dependence of t has been neglected from this result and. 
as such underpredicts the dissipation for large Mach numbers. 

The sensitivity of turbulent viscosity model to the selection 

of C .is the final parameter examined in this studv. Figure 27 
turo ^ 

shews that an increase in C . gives an increase In docav rate. This 

turb 

paiametor has been treated as having secondary importance since only 
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qualitative results can at this tine be predicted. More reliable 
turbulence data is necessary to assure the proper nodel for the tur- 
bulent viscosity or a proper value of C . . 
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Section VI 
CONCLUSIONS 

A theoretical study of the stability of flows within coabustlon 
chambers with evenly spaced baffle configurations was presented. For the 
first time, a stabilizing influence for baffles has been properly pre- 
dicted in an analytical model which incorporates the influences of a 
concentrated combustion source at the injector, a "short" nozzle ter- 
minating the chamber, and mechanical energy dissipation at the baffle 
blade tips. Two separate combustion chamber geometries are examined 
with this model and produce the following results: 

1. The addition of a baffle to a combustor in many situations 
will improve the stability of a chamber. 

2. A fluid dynamic loss created by the effects of viscosity 
and turbulence produce tne damping mechanism of the baffle. 

This energy dissipation occurs locally at the baffle blade 
tips. 

3. Without the effects of mechanical energy dissipation, wave 
alteration produced by the addition of a baffle to a com- 
bustor cause a destabilizing influence. 

4. The baffle is most effective in damping the traveling 
transverse modes of oscillation. 

5. 'n improvement in combustor stability is generally achieved 
with an increase in baffle blade length. 

6. Longer baffles may be required for combustors which contain 


an Increased mean flow. 



The addition of a baffle to a combustion chamber depresses 
the oscillation frequency. 

Baffle configuration has a secondary Influence on the 
stability of the chamber. 
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Figure 1. Several baffle configurations and baffle blade shapes. 
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Figure 2. Aerojet-General experlnental data of decay in decibola/eycle ve. 
baffle bladu length. 
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Figure 4. The two dimensional rectangular chamber. 
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Figure 5. Axial velocity and velocity potential at th« main chamber-baffle 
compartment Interface in a chamber with a two compartment baffle 
with blade length 
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FlC'ure 6. Axial velociry and velocity potenCia] at the 
inain charibef-baf f le csaspartnenc interface in 
a chamber with a three compartment baffle 
with blade length Zg*0.3. 
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Figare 7. Polar coordinate syscaa at baffle blade tip. 
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Chamber conditions: 


L=1.5; K-=0.0; 



Figure 10. A streamline plot of the unsteady flow within 
a baffled combustor (first transverse mode 
in the main chamber). 
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Figure 11. The geometry of the three dimensional baffled 
chamber . 






Figure 12. Axial velocity and velocity potential at the main 

chamber-baffle compartment Interface 'and r=l.), 

in a cylindrical chamber with a three compartment 

baffle of blade length z *0.3. 
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Figure 13. The comparison of Che solutions of axial velocity at the main 

chamber-baffle compartment'^ Interface using the series expansions 
summed without Cesaro summation. 
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Figure 14. The comparison of the solutions of axial velocity at the main 
chamber-baffle compartments Interface using Cesaro summation 
of the series expansions. 
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Figure 16. 


Prediction of the frequency depresnion vh. baffle biade lentth in a 
Two dimensional chamber. 
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F:<.gure 17. Mean flow Mach number corractlcna to the frequency ve. baffle blade 
length. 
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Figtire 18. Decay in declbela/cycla va. baffle blade length in a chamber which 
haa no combustion, nosnle or mean flow effects for various wave 
amplitudes . 
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Figure 19. M*ian flow, combuotion a.id nozr.ltt influence# on the 
prediction of decay in dcclbele/cycle ve. baffle 
blade length. 
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FlRure iO The effect of number of baffle conpartmenta on the damping 
in decibele/cycle va. baffle blade length. 
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Figure 21. The predictions of damping In decibeln/cycle va. baffle 
blade length for various oscillation modes in the main 
chamber . 
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Figiirt" 2 ?. Thu decay rate ( X ) vs. baffla blade lent^tii for the flret and 
third transverse mode In the main chamber. 




Figure 


3. Normalized frev^visncy ft. baffle blade length In a three 
vilmenslonal chanber. 
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Figure 24. The effect of wave amplitude on the prediction of decay 
In declbela/cycle vs. baffle blade leng> 
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Figure 25. The standing wave and traveling wave predictions 

of decay in decibels/cycle vs. baffle blade length. 
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Figure 26. Mean flow Mech number corrections on decay 
in decibeis/cycle vs baffle blade length 
(wave amplitude *‘0.1). 
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Figure 27. The effect of the turbulence coefficient on the stability 

predictions. 
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APPENDIX A 


The set of constants which appear in Equation 7, 8, 19 and 20 are 
defined by the partial differential equation and the appropriate 
boundary conditions. For the two dimensional chamber these constants 
are defined as: 

Baffle compartment - Two Dimensional Cha.aber 


u)M + 1- (M^ - 1) ((ntirN)^ - co^) 


2.B 


OJM + (H^ - 1) (faTTN)^ -u>^) 

(1 - M*) 


S 


•I- 


Bi „ + M [yad - e"^”^) - i] (OJ + B^ ) 


.B 


Bj g + M [ynd - e“^“") - l] (u + Bj^^) 


] 


Haln Chamber - Two Dimensional Combustor 


oiM -t- (M^ - 1) ((mir)^ - O)^) 

i.c ~ ^ ^ 


d - MO 


oj M - + (M^ - 1) ((mrr)^ - 

d - mO 


C = - 
c 


B, + M(y - 1) (u + B, )/2 


Bj ^ + M(y - i) ( + J/2 J 

,C 9 •. 


For the three dimensional cylindrical chamber these constants are 
represented as: 


PRECEDING PAGE NOT FHifE0 



A-2 


Baffle compartment - Three Dimensional chcF^ ~- 




Maln 


.c 


B 


2 


jC 


C 

c 


oM + 

+ 


- 1 ) 

<■ 
' im 

- 0)2) 


(i - 

M^) 




M - 

+ 


- 1 ) 


- w* 


(1 - 

M^) 





g + '.1 tyod 

-iuTT. 

- e ) - 

1 ] (u + Bj g) 

» 

g + M [ynd 

-itdT. 

- e ) - 

l] (o) + Bj g) 

iiamber 

- Three 

Dl'.ensional Combust v*: 

(»jM + 

^ + 

(M-- 

- 1 ) 


- u)2) 


(1 - 

M*) 




UM - 

+ 

(M^ 

- 1 ) 


- U)2) 


(-• - 

m") 




r "s 

^ + M(y - 

- 1 ) 

(u + 

«^c 

)/2 -j 

Lb,. 

+ m(y - 

- 1 ) 

(U) + 

»2.e 

)/2 J 
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COMPUTER PROGRAMS 

The coaputer programs usc<i to calculate the frequency (<^) and the 
decay rate (>.} ^‘f oscillaclcas for baffled cosbustors aodeled with the 
iofluences of gaia/i-:i>s en<I~vall boundary conditions, nean flow, and 
baffle Lip energy dissepation are given in this section. These prOc^^raas 
ace codec In Fortran IV using coanlcx arithenatic. Two pr> graas are 
written fet' the two-diaensional c!<aabri (FiS>B) and the three diaienslonal 
cylindricul chaaber (BAFFLER. Terainating each prograa is a saaple output. 

Two Diaeasional Chaaber - Coaputer Prograas TWDB 

The input paras, eters to this prograa are devlnei asr 

1. MC - Maximum nuaber of series terms used tc represent the solu- 
txon in the aain chaaber. 

2. M3 - Maximum nuaber of series terns used to represent the colu- 
tlon in the baffle coaparcaents. 

3. ALENGTK - Nondiaensional chamber length (L*/R*). 

■». ZB - Kondiaensional baffle blade length (2 */R*) . 

O 

5. T - Noudimensional baffle blade thickness (T*/R*'. 

6. MUB - Total nnaber of evenly spaced baffle compartments. 

7. GAMMA — Ratio of specific ^ats. 

8. ErSllON - Oscillaticin vave anpl 'tud-.-. 

9. *J<ACU - Steady state nean ilow Mac:, number. 

10. - The iransveise t&'de numb.^i cf tlit- Honir.ating term in 
the nai.T chamber solution. 

11. 1DM.V .Maximum number . » i'-i-f-.tion.s of tiie successive apprexi- 
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Ocher initialized paraaeCers (not read in Che progm as data): 

1. Vise*- Mondiaens local molecular viscosity (U^*/ p*a*R*). 

2. CTURB - Coefficient in Spalding's effective turbulent viscosfty 
L- /del. 

:5. AS - Interaction index paraaeter of C.occo's ciae^lag aodel 
unsteady combustion 

4. TAU - Time lag of Crocco's tiae lag aodel. 

The input data for progyaa TWDB is read on three cards and has the 
following format.: 



Columns 

Variable 

Tvge 

Card 1 

1-5 

MC 

lOTEGER 


6-10 

MB 

IIJTEG^ 


11-15 

lOMAX 

INTEGER 


16-20 

MHAT 

INTEGER 

Card 2 

1-10 » 

ALEHGTH 

•jBCIhAL 


11-20 

AMACH 

DECIHAL 


21-30 

GA.M«A 

jECIMAI- 


31-40 

EPSILON 

DECIMAL 

Card 3 

1-5 

HUB 

lyiELER 


6-15 

ZB 

DECIMAL 


■ 16-25 

T 

DECIMAL 


« 

A discing of the program ar^ a ‘sample output is fontained on the 


iollcving pages. 
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Three Diaensioaal Cyllodrical Chaaber-Progran BAFFLE 

The inpuc parameters to this program are def -ned as: 

1. MC - Maximum number of Fourier series turns in the main 
cnamber solution. 

2. LC - Maximum number of Bessel terms in the main chai^er solution. 

1. MB - Maximum c-rmber of Fourier terms in the baffle compartment 
soutioa. 

4. LB - Maximum number of Bessel terms in the baffle compat*tment 
solution. 

.5. ALEMGTH - Combustor length to radius ratio (L*/R*). 

6. ab - Hondimensional baffle blade length (Z3*/R*) . 

7. T - ^ondisensional baffle blade thickness (T*/3l*). 

8. MUC - Number cf evenly spaced baffle cavities. 

9. GAMMA - Ratio of specific heats. 

10. EPSILON - Oscillation wave amplitude 

11. AMaCH - Steady state mean flow Mach number. 

12. MHAT - Dxainating transverse node of the oscilLaticn in rhe 
33 in chamber. 

13. LHAT - Dooinatiag radial mode of the oscillaticn in the riain 
cnamh.;r . 

14. IDMAX - Maximum nuab.r of iCtrations of the successive approxi- 
mations of the solution. 

Other initialized parameters which are not Input as data: 

• 1. Vise - Moleculcr viscosit r (u */.'*a*R*) 

v 

2. CTURB - XTcef f icienr ic Spalding's effective turruient viscosity 


model . 
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3. AK - Interaction index. 

4. TAU - Tiaie lag 

3. MX - Variable which detersiaes the type of oscillation to be 
considered. MX = 0 will specify a standing wave solution and 
MX = 1 will yield a traveling wave solution. 

The input data for prograo BAFFLE is read on three data cards: 



Coluans 

Variable 

IZEl 

Card 1 

i-5 

MC 

INTEGER 


6-10 

LC 

t:<teger 


11-15 

K£ 

INTEGER 


16-20 

15 

INTEGER 


21-25 

ID.MAX 

INTEGER 


26-30 

MHAT 

INTEGER 


31-2j 

LHAT 

INTEGER 

Card 2- 

0 

1 

ALEliGia 

UfCIMAL 


li-20 

AMAC-i 

?ECr«AL 


21-30 

GAMMA 

DECLMAL 


31-40 

EPSILOX 

DECLMAL 

Card 3 

1-5 

MU 8 

INTEGER 


6-13 

28 

DFC IMAL 


16-25 

T 

DECIX.AL 


A listing and a sample output is c 'ntaired in the foliowiai; pages 




PROGR*** T«Ol»( lhPUT.OUTPOT.TfcP£S=lM»UT.T*Pc6=OUlPUT) 

COKPtO »USKtA8t«*d3»Cn K83*SURv,C{«BMUt 

lALAMK.C6S«C^l.C8afCCBIttC82*CCC»9SV>CCS*BBl »Sd2*CBl »EB2»6CC-»Gl tOt 
2OH€Gl*At!Nj«UN0ZtBA:>E*SUMZ«SUMUSUH2«VI.2«VL3»AHu*A tOMEGteK 
OlMENSION TTEi30- 3C«lt)l »AKUf 30 1 101 «m< 1 1 ) *Zf 1 1 1 *A< 101 *0<«i 

0ATAOi<Jl.i = ltlU/.e212'#101B375SA..057&166583112A..0934003V82782S. 

1.12052Cl69GH»32».13fc^2*«9o09S628».lAU9J706*5287S*,136«i2AV0B^628. 

2. 120520 l69Gl«32».093‘.0039827825t.OS7616feSe?112A#,02129lo!837b3*/ 
OATAlZ<n« I = l.lil/.00 79573l99S2S8..0-V6VlC07703067,,12291663fl«.se. 
1. £30765 3AA9«*716«.36018<'793<^1911«.S00000090r 0000*. &39S1S206S5U89. 

2. 7692i^S6S53S26A,.8/ •083363285*2., 9530899229693 J.. 99, ?0*2600A 7*2/ 
fNHUlS IHCtOOE NO. Of SEHIES T£kMS-»*C AND H8»COM0uSlOR 
nONOIkEKSIONAE LENGIh-ALEHGTH, baffle blade LEn6Tm-ZB» 

c**»* bmeele blade rmcANEss-i.NO. or baffle corpabtmekts-hub* 

£•••« ratio Cf SPECIFIC MEATS-GAMMA.rfAvC AHPU 1 lUUE-EPSItON* 

kEAN FLOta 8ACH NO.- A8 ACh,TANSVEN 5E PCDE C^APACTE^i-MMAt . 

£•••• and fASlKUM NO. Of HEHATIONS OF SUCCESSIVE APP«OX. Of SERIES 
(iEAUIS*l*0) MC.*‘e»IOnAX,«hA.T 

1*6 FORPAT^AIS) 

R£ AO ( 5 1 i * I : ALENG TN . AH AC<n. GAMMA .EPSILON 

l*ti fOH:!A. t*r lO.SJ 

k£A0!3.>*?inu3,ZB«T 

1*1 rCkMAT ( IS. 2f 10.51 
aC=XB=Z8 

C***« these FaRAHETERS are for the DlSSIPATf IN CALC. 

CTL*Rc=C.03* 

VlSC=1.0E-05 
OiSfV=5C*T tviSC) 

0iSFT=£0«T (STUR8I 

^•e«» uHSTEACt CO^^BUSTION iNPuT PAHAHElEkS N AND TAU *HE DEFINED 
AN=(CAFKA«!.)/*./GAHMA 
1AU=1. 

WR1TE(E.*00> 
ttHlTEfE.AOOi 
«r(l T£ (E .*0 1 • 

*0! FOkHaT(«0ThE following CALC. ARE MADF: FOR A 2*0 ChAHBEk«1 
•R1TE(E.*00> 
aRITE (E.*00> 

*00 fohmah»o«. 5o«ih*»» 

WHnE(E.*02> 

*0< FORHAr(*OTHt CnAmBEH ha^ The FOLLOhIvG CONFIG.*) 

HRiTE(E«*03> AlENGTh.AHaCh.GAMKA 

a 03 fOKHAT (•0L£NGTH=*,F 10.S«3X.« hEAn FlO* HACH NO.=«.Flo.St 
13X« / •• RATIO OF SRtClFIC H£AfS=-».» 10.5) 

WRITE (E.-O*) h,,i6.ZB*T 

*0* FORRAlMClHt baffle consists OF *.I5.* COhPARTMEmTS wITh LE»WGTh». 
lFlC.5»/«‘ AND BLADE rhlCi‘NESS*»f 1 0 .5) 

*r1TE(E.*05> hhat. epsilon 

*05 FORMAT (*0T«t OSCILLATION CONSIST OF A COM INA T ING* . I 3« • TRANSVERSE 
1 NODE**/.* and IS OF WAvt AHPLlTvFC-^.Gai.l*) 
wR 1TE(E.337I Vise 

337 FORMAT (•5MOLECULAR VlSCOSITVr*.G?l . 1*1 
WRITE <E. ‘ 071 CTL'RS 

*07 FORMAT ••OEFFECT IVt TuRBULcNT VISC COEF .s«,021 . !*• 

WRITE (E.*08> 

*08 FORMATi*CTHE FREOuENCV ITERATIONS AFE AS FOLLOaS*) 

UU 100 ML=1*MUB 

OPsMG 

OR2:PUE 

w;»Cfi/CR2 

V2slOR-l.)/OR2 


AL Page is 

PRECED2^G PAGE BIANX M03C IfUaSl * (WJALITE 



00 190 HSTAH=i,H 8 
00 100 M=l*MC 
B1=3.1A159265359*HMSTAR-1 
B0-dl*El 

lF(e0.tT.1.0t-06l GOTO 101 
r*2=3.l4il5926S359» ( (MSTAH-l . • »HU8* (M>1 . n 
8AS1=(S1»I(B1*V1I-SIMB1*V2) l/t2.«bl) 
8AS2=»5lNlB2»Vll-SI*«(B2*V2n/<2.*B2l 
TTE IMST AR,M*MU1 -<3AS1 *8AS2 
GOTO 100 
1m1 EP=1. 

IF(H.EC.l) EP=2« 

TTE(»«STA«,«.MuI=EP/*2.*MUB» 

100 CONTIM.E 
999 COHIINLE 
S60 CONTINCE 

€•••• first -c-uess of frequency 

ONE 3-Mt- i T • 1 3 . 1 1 59365 359-2 . • ' 3 1 
00 300 N=i.NC 
300 bMulN)=CHPLXf0.0«0.0> 

RHuCHhAT*l)=C» PLX«> «6. v.Cl 
*'i-3.!«c1592633S« 

CIsCPPLXiO. 0*1.0) 

MSsMhAT •! 

0(l)=C)'PLX ( .*,0.0) 

0;2):CNPLX(0*0«1.0) 

0i3)sCPPLx I..O.OI 
O(n)=C?*PLXI0.C«-1.9) 

10=0 

09 COnTINCE 
uEl =0 • 1 

SUHX=0’PLX (O.OtO.O) 

AN' aN*n,-CExP(-Cl»OMEG« TAO» 1 - 1 ./GA»<m» 

ALANl=-GAMNA»Cl«ONEG«A^.AC«*ANN/n .•t-ANMA*AMACH*AMXCN*ANNI 
£PH=l* 

JFTNf-AT.EO.O) 

IFdJ.EO.Ol GOTO 500 
EPS=2. 

00 502 NST = 1.’ 5 
0«=C**ST-i. • "S'NO-l 

SSK=CNEG*ONtG*AN*CM*AMACH. | AKACH* AKAC". ' 1 , ) • iOR’'OP-ONE5»OHEG) 
BSKsCSCRT ;SSAi 

Abl = l0t'EO*AMACh«eS»<5 ' ( 1 ,-AKACh*AHAC»- » 
Ad3=(0>'Efe*ANACH-aSK)/U.-ANACH»AMACH) 

CMs-TCMABI-ALANII / CCI*Ae?-AtANl I 
E8!-CExP(CI*Abl*2e) 

Eb2=CExP(C»**B2*Z«» 

EB3=CI* f Abl •101 •C*»*A82*EB2 )/(£«» *CN«Eb2) 

SUMYsCnPEX <6.0.0- 0) 

00 503 H=1,NC 

IT (N.EC.NS) GOTO 503 

TUHzQ.O 

00 SON Nu=l»»*U0 

50-<. TUMsTU^^TTE (NSTfMSiNGI^TTE (NST.«.»*U» 

SUNY»StNY*BNU(«» »TUN 

503 CCNTINLE 

SMMX=SCNX*A.«M'ja*E83*SUNY/tPH/EPS 
502 i.PS-1. 

5u9 CONTINlE 
22I- IN=0 
215 CONTINtE 



ANN= AK« ( i . -Cf XP < -C 1 •out G« T Aun - 1 . /GAMMA 

AI-AMI =-C-AMMA*C I *CME6*AMACH*AAiN/ 1 1 . ♦6AmmA«AKACH* AKACh«AMN» 
A8 sAMACh*cGAMMA-! .!//. 

ALAMK=-Cl«OyE&*A«/ ( | « •A mACH«A 6 i 

CbS = CSC«T (CMEG*OMt&»AHACM*AMACH» I AHAC»i«AMACM-l . I * <KhAT«MMAT*Pl«PI- 

l0ME0*0“h;,) » 

CblsJCMfcD^AMACM^Cti? »/••! -AMACh-AmAC>-J 
CB2= (C»’E.&‘»/.*<ACH-''SS* / ? i .-AMACi1*A#<ACM 
£C8l=CCX?(CI*CSl*JZU-ALtHGT»-» I 
PCBZ=Ct*P«Ci*CH?*«fd-Al-tt<:&TH} = 

CCC-- (C i*Cfil-Ai A^r j/jci*Ci^Z-al*^''« 

BSV=Ci* (C8l*ECiJl‘CdZ»CCC«FCr-2)/(tCu:*CCC*£C0Z» 

SU«Y=C^‘Pl X 10. 0*0. 01 
tPS=2. 

00 5S^ M=l,MR 
OP=(M-l.J*PI*Mob 

CCS=CSCMT iOmEG^CmEGoamaCH^AwACH. i A^•ACH• am AC«- 1 - » • ;oh»or-omeg*omeg> 

u 

681 = iCMILGyAMACH^CCS) / 1 1 .-amach«a»:aCm) 

0BZ= ;0*-E.G»AMACM-CC3t / (1 .-am<C«*'Amach» 

eHl=C£:xP(CS*88l*ZriJ 

£02=CExP|CI*BeZ»28> 

8CC=~ it I *BB I - AL am : ) / ; C I *882 - AL AM I « 
xgH-O.C 

00 j51 hu=l5MU6 
551 TUH=YV;»’>lTTt«M,MS«Mijn**2 

SUMY=SLMr.<..«Hue*TUM*CI* (861*£81 <ea2»eCC*EB?i/£PS/tPH/CtBl»6CC*t83 

}} 

550 EPb=i. 

Gi =HSV- StiK*^-SuMx 
6=CABStGl ) 

IF(lN.NF.O> GOTO 2P*. 

223 OMEGl=CH£G 
m£F=C 

GOTO 225 

22A lF(6.Lf.*FFI GOTO 223 
lFilK.LT.Ai GOTO 228 
OMFCsC^EGl 
0EL=r>£L/2. 

Ir <0EL.LT.i.F.-03) GOTO 310 
IN=1 

GPT3 225 
22» in=;k»i 

22= 0«tb:=0>'t.GJ ‘Ui !K)*0EL 
GOTO 215 
310 CCNIIME 
10=iC*l. 

•KITE ,'C.Al£> ir.C*^Eb 

Alo FpK5iA! l«OlTtPATlGN ‘.U.* COmPlEx FKEO.s* . 2 G 21 • Ia» 

(IC.CT . IDMAx) goto 900 

C«A«*> IT2RAT10K C^ MAlK CMAMjtM 5tBIES COEF. 

ALlNJ = ALAM 
AINCZ=ALAMN 
EPP=Z. 

00 200 MP=I.**c 

IF 1MF.E0.M5J GOTO 200 

0K=(MP-1.)*PI 

b*St=0J'E0*CMEG<>4MACn»AMAC*'» ( amACh^amach-i . ) • (0«*Oft-OMEG«OMEG » 
3ASE=CS0RT lOASE) 

CBl=(OMEG»AMACM*aASt)/( j.-A“ACH»*MACHI 


I : 





€H2a <0»>c.0* AM«CH-b*S^. ) / ( 1 .-AHACH«A»»ACMI 

CCC= ‘(Cl -Cei - AI.N02 ) / « C I •f.B2-ALN02 1 

ECBi=C£xP(Cl*CEl*(2U‘At.E.N6rH}l 

CCe2-'CexP(CI*CB2*(2tl>Ai.ENGTH)> 

C«S=CI«JC81*tC81»CB2*CCC«ece2»/(tCB!*CCC«£C3»l 

SUHY-CPPLXCO.O.G.O) 

suHx=0’PLX(0.o««;.o) 

tPI-E. 

CO iol HST=1»M8 

or=(mgi-i.;*pi*hu0 

base -C£ OPT (OMEG*nMEb«AHACH*AMACH*(tH«C^<'AMACH>|. >*(OR*OK-OMEO*Of4E6 

1»» 

B81=(0PEG*AM*Ch»S*Stl/{ I .-AMACH«A mxch j 
B82=l0P£6»AM4CM-6AStl/Cl,-4KiCM«ftHXCM» 

BCC=- (Cl*ti 0 l-ALlNJ|/cCI«a«i 2 -ALINji 
EBl=CEXPlCI*B8l*ZB) 

EBif=CEXP(Cl*862*2a» 

TUHsO.O 

DO 202 HU=1*MU6 

202 TUM=lui»«<TTtCHST.MPfhU> J**2 
BAS£-Cl*(B3l«E81«B82*8CC*ER2>/(EBI*eCC*E82) 
SU»<X=SLMX*<i«*Pua»BASE*TU«/£PS/6PP/CeS 
SUM2=CPPLXJ0. 0.0.0) 

DO 203 M=1,MC 

IFlH.EC.hP) GOTO 203 

TUHsO.O ■ . 

CO 20A MU=!.Mua * 

20A 1UM=TDP»TT£IKST.KP«HU;'*TTE (MSI .MtMU) 

SUH2=SCM2 *HMU (M) *TUM 

203 CONliNLE 

SUMY = SEMY*SUM2*<i,*MUb*8ASF/£PS/£PV/Cas 
201 £PS=1. 

bMU (HP » 3 SUMY/ ( I .-SVMX) 

200 f.PP=l, 

GOTO 8S 

900 CONTINUE 
bWITE (C.91S» 

VJ5 fOKMAI (•OT»it MATPIX Uf T«E M*|N CbAMflFP FOURIER COEF, = *> 

00 9H Msj.HC 
uRITE (« 

910 FOMHAT(« •,13.» •.2021.1a) 

911 COlyllNLE 

00 <<01 M«j=l.Mua 
00 901 Msl.MH 
SUHX=C»'P>.X(0. 0.0.0) 

(JO 902 HPsl.MC 

902 SUMX=SV.**X.bMOtMP) *T»t (M.Mt’.MU) 
tPsl. 

IF(M.Et.l) tP=2. 

BASEX3FP/(2.*My8) 

0M=Mb-P.l. 

OKI =KB 

SUMX = SLMX*OM(rOBl 

901 AMU(* .M;J =i-<1*/8/SEX 
DO 913 MU=:.MUn 
•nlTl ( c.91a) My 

CO 912 Mrl.Mii 

912 ■K1TE<E .RiOfM.AMyiM.Md 
91j COnTINcE 

9IA format (•0)«t MATRIX OF FCmRIEP COEF. FOR Tm£ BAFFlE CO»*P.».l3) 

9er COM SUE 



MUBX=MLE~l 

00 778 MU=l*Mi;8X 

SUMl =SLH2 =CHPlX O .0 .0 

YlsMU/PUa-T 

Y2=MU/t'UB»T 

00 779 Hsi.MB 

VL1=(M-1.i*F1*MUP 

VL2=*HC (M«HUI«C0S(VL1*Y1I 

VL3=*M- lH,Mb*n«C0S(VLl*Y2» 

SLMl=SLHl*Vi.2 
77S SUH2=SCP2*VL3 

*<MU) = (SbM2-SuMl)/SUHT t2.*Tl 
778 CONTINLE 

£•••• dissipation '.ALCuLATION 

0M=PEAL (OkEO) 

AHULT=CAHMA«S0HT (OM /2*>/2. 

SUHsO.O 

7£=ACOS(2B/(2B*T)) 

00 339 Pb=l*NUBX 
ABX=CAES(A(NU) ) 

A8X=ABX«ABX 

B*EPS IL0N*EPS IlON* A8Xv'T/4 . 

SUH1=0*0 

00 336 1=1*11 

tTA=IE/2. tZU J*CPI-TEI/2. 

Vl=SiN(ETA) 

338 SUMlsSLMl* <0ISEV*0ISFT« U AMACH»AMACP»B*Vl*yl » I)/V1 

S 0 Ml=AEx«SUMl*(PI-TE»/ 2 . 

339 SOH=SbP*SUMl 
SOH=AHtLI«SUM 
OHC=HEA). (OMEGI 

£•*•* oecay pate calculation 

B1NI=0*0 

SJ“tl=0,0 

tPK=2* 

DO 821 H=1 *Ni 
0B=PM<N-l.-‘*KOB 

bSKsCSCfO (OhEG*C‘'EG*AMACn»AMACH* , J* «0«*0P“0 pc0»0»'EG» 

1.' 

ABl = <C^EG«' AM4C‘.»BS*< ) - ( 1 . -amach«AM*Ck > 

A82= «0»<E&*AM4rH-3SK)/ ( i .>AKACH*’ i«ACP» 

CM=-fCl*ABl-4LAPJ )/ (C1*AB?-ALAH1I 
BASE=C£ xPCi'.' >11*26* •CK«CExPCCi«Afci2«Z8» 
dAS=CAES(EASE) 

BAS-HAS* jAS 

«81= ( 1 . 'Ch) •L.*0»«EG*APACVi*CI* I AtJl *Ae2*CH» 

BASM=C^6S (faol I 
8A_« = eASf<*3ASW*C-A»*'lA*G4>"*A 
SU"m1=C .0 
DO 622 1=1.11 
ETA = 2d*2(U 

8SA=Cf P(ri*ABl*F TAJ *Ch*CExP|C1*A«2*LTa) 

BSi=CttS<BSi<l 

BH?=CI*I i'l'CEr.PtCI 'AB1*ETAJ .C*'*Ar2*Lc*P«CI*A^2*cTA: ; 

BS' =C A^«A* (CI*ONEG*«a AMACh*082 » 

BS =CAt<;<t»SM 

BSsiJS'^S/f » /G*»4f*A»GAhMA*0>*G«0M »ijSl*BS»/2. 
t ?2 SU*^*1 =£uhk. »dS*»< P*<fc/t*-? 

»4H=j.O 

0 823 Nt=l*Mui 

HAHl =CAbS(AKU(M.MU> > 



V 



623 RAM=RA»< »RAH1«RAM1 

BlNl=BlKl«RAH«BASR*tRH/BAS/FLOAI iMUB>/2. 

SUHY 1 =£UMY 1 ♦RAH*£PH«SOM»t 1 /FLOAT <MU8 ) /2. 

621 EPMsl. 

6NZ1=0.0 

SOMOUC.O 

EPM=2. 

00 8A1 Hsl«NC 
0R=PI«(M>1.) 

8K=CSGRT (OHtG«OPEC«AMACH*AMACH*|AMACH»AMACH-I.I*iOR*OR-OHEO*OHEGJ ) 
CB1= (0PEG*AMACH»8K» / ( 1 .-A«ACM*AMACh> 

C82= (OPCG*Amach-BK » / ( 1 .-AMACM»AM ^Chl 
C‘«=-lCl*Cei-ALAMM)/»CI«CB2-ALANNj 

BASE=CExP(Cl*C81*(Za-ALENGTMH*CM*CExP<CI*CB2*(23-ALENGIM» » 

BAS=CAeS(PASE) 

BA5=6A3«BAS 

8BI=(1.*CM)*CI*0mEG*AMACH»CI»(C8I*CM*CB2* 

BASR=CABSt83n 

basr=6asr«basr*gapsa«gahha 

SUHVlsO.O 

00 BA2 Isltll 

ETA=ZB* ULENGTH-Z8) » I » 

BSK-CEXP(CI*CBl»(ETA-ALtNGTll l•tH*CEXP|CI*CB^•(£TA-ALENOTH) ) 

8Sl=CAeS(BSR> 

8B2=C!«{CBI*CExP(CI*C81*fETA-ALEN6TH>»»CM«CB2*CExP(Cl»CtJ2*( 
IETA-ALENGTh»> ) 

BS^=CA►MA•(C1•0»*EG•BSK*A^*ACH•BB^I 
ANZB = A»' ACH« ( 1 .«C-AHHA) /2./GAMMA 
BS =CAES(BSK) 

BS=8S*ES/2./Gamma*GAMma*JBG»Gh6«6SI*BS1/2. 

J»<*2 SUMVl=SUhvi ‘BS** ( 1 1 • ( ALf NGTH-Zbl/eAS 
Rahi=CA8S(BHuiM} ) 

tiNZ I =8A Z I *R AM J AM 1 *b A?,H*EPM/&AS/2 . 

SUMUl=£bM01 *RAM1*RAM1*'SUMV1*EPM/2, 

8<il LPM=1, 

iilNlR:ElNi/2. 

riNZlK=eNZI/2. 

kn=ak* o .'- cos(Ohg*taui »*amach 
•j!N!RrrN‘*bIBIR 
BKZ'iP=ANZe«BNZIR 
SUKlN=6KZlR-BlNlP*'iLiM 
«MITE(E*30) BINIR>BNZIR«SUM 

30 FO^HAT «»0C0M8. IN JuT*tGl3.6.« NOZZLE E*7 . =*.G1 1 . 6. • 0!SS! 3.6 
I i 

SUHY1=S0MU1 ♦SUMVI 
OECYWIsSLMIN/SO’ .1 
»iMlTE<E.12S» OECtWT 
i:S fOHMAT (•OOECAt R A TE=« .G21 . ! A > 

OECTMseA.s^^-^CECrPT/OMeo 
KMlT€(e*13Ci OECTd 

>:0 format (*OOECAy in OLCIHELS /CyCLE = *.G2J .1<,» 

END 



THE FOLLOmISG c*lc. *ke m»oe foh a 2-C chamber 


THE chamber has the FCLLOhING COHFIG. 


tEH&Ihr l.SOOOO HE*n 
RATIO CF SPECIFIC HEAIS= 


I LOW MACm nO.= .10000 
I >20000 


THE baffle COKSISTS OF 
AMO BLADE THICkmESS 


2 COMPAMTMEMTS WITH LENGTH 
.05000 


.SOOOO 


the OSCILLATICM consist of a DOMINATING 
AND IS OF WAVE AMPLITlDE .10000000000000 


transverse mode 


MOLECULAR V'SCOSiTT= 1 .OOOOOOOOOOOOOOE-OS 

TUNBulENT Use COEF.s 3.AOOOOOOOOOOOOOE-02 
THE FREQUENCY ITERATICNS AME AS FOLi.OwS 


ITERATION 
iteration 
iteration 
iteration 
ITERAT ion 
iteration 


1 COMPLEX FrEO.= 

2 COMPLEX f PEC. I 

3 COMPLEX .'RtO.r 
<• complex FNEt.= 
5 CCMPlEx fneg.= 
o COMPLEX FRCO.= 


2.21SO3O1S3S9O0 

2.1 3221 76S3S899 

2.1 3221 76S3S899 

2.1 3221 76535899 
2.13221 .oC?5899 

2.1322176535899 


THE MATRIX OF ThE MAlh C-iANBER FOURIER COEF.s 

1 ”<•. 18151 85260 01 t2E-13“6. 780 na669 7729S9E'*13 

2 l.OOOOOOOOOOOOC 9 . 

3 ■! .763S3I965O8053E-1 3 2.07«»276l6893762C“15 

A U3S9<iS607o 026 . *33 1 9N Jd637603E-C3 

5 1.5*2157979509<.2£~13 2.*0557I2*57«.6S6E-1S 

6 6.!.7652700865522t-02 9 . 81 393901 68* 786E -0* 

7 -* .6378e*506**2Ut.-l* I.S651 720*97537*£--16 

6 '3.3**2781 ■ u- 1 37*1 -o^-S . **3*26U263v»55F-0* 

/ 1 >0*9699766o066et-l 3 1 . 1899S8o53o9389E~15 

10 2.60*81 j010J6568E-U2 3.52 78** j799J>i21E“0* 

11 1 •03*999£*.23912eE-l* 2.6?3596U 995957E~1 7 

12 -1 .90870 17 55 1683 IE -02-2. 50525* 75556*605-0* 

13 B.6091980S*70*82E-1* 7 .<,6803**9 76 1 95 6 

1* 1 .*7J221*/SJ0077E-02 1 . 88906*91 0*83 78E-0* 

15 2.3«8510t05J25*3E-l* I .290093625 1 5302E-1 6 

1 6 -1 '180117tl6299*3E-02-l.*a6138'.5B2d587E-0* 

17 e.OI 1833C9931*82E-1* 5.39e7902.>036606E-i6 

18 9. 72n6U272a7eOE -03 I .206630' 679* J59E-0* 

19 3.2666*9*08o8263E-l*-2.8i 729502*5e 1 09E-] 7 
2 u - 8 . iR *6 7.e50HdS3ec - 03 -) , 00 39*02ol f>-ju 6 7E-0* 

21 7.6*S75je2651355E-l* * . 0993 72 y 5356557 E -1 6 

22 7.013*7203159*331-03 t. Si 780c J 6 ) 922*9E -05 

23 3 . 69*856 8** 0660 * E - I *-* . 2900 7***9623 73E-1 7 
2* -6.09790 7** 6 76 790 1 -0 3- 7, 3* 3<»9 75055932 IE -05 

25 6.96009061366566E-1* 3. 1 > 3*2732556S93E- J 6 

26 5.J67257Ei-.92653E-03 6 .* S 65*C63C96566E-0S 

27 3.801 y*57 !222J83E-1*-5.37*?J9COJ3*70S£-17 
26 -*.7 7* 10 69 0661 1 90E - 03 - 5 . 67 1 05 7 1527*1 19E-05 

29 6.687*330/1 7''365E-1* 2 .60 1 59**926669?C-1 6 

30 *.28562656*09*73c-03 5. 062 1 02932891 SOC'OS 


'.lOOOOOOOOOOQOO 
-9 . 667*9999999998E-C2 
“9.5312*99999 Vv95E-02 
“9.5312*999999995E-02 
-=.5312*999999995E-02 
“9.5312*9999V9995E-02 







THE MATRIX OF FOUHIER COEF . FO« IHt BAFFLE COMP, 1 

1 .68‘<430007b3Bt9 ».6239I8M79«»bS5E-0A 

2 .266632IS10A316 -2.33152099A2ASS3E-03 

3 “e.82S65Te3269'>A9E-02-3.'‘eMfl09«<.7lao«AE-OA 
A A.686I26393I32SSE-02 2.99882 i29|60«i0AE-0A 

5 -2.982Ca912A8»28AC-02-2.3090038A9A<f334E-0A , 

6 2*093Ar03805At92£-02 I .792A88A902372AE-0A 

7 >I.S6393908730201E-O2-i«4I0ASA2<.2eo<.38E‘0A 

8 1.?201i)1A;v<>1971E>02 1.IAA179L8163099E-0A 

9 -9.832lS20697b0«*6E-03-9.€,8eA6217251t>26E-05 

10 8.13O1291820O8O3E>O3 8.02320923332bA9E-05 

11 -6-.87193021073A5eE-03-6.9066260e-*58939E-05 

12 5.9263351 30989ASE-03 6.0515127 V 19>485E-C5 

13 -5. 220060993 737 12F -03-5. <.076d933505575E-0S 

19 9.7222e7S55bl305E-03 9.96181907292996E-05 

15 -9. 97718760309159E-03-9. 779762199 16039E-05 

16 S.07390708379687E-03 5.5621027525J829E-05 

17 -3.19S126792U673E-03-3.9199692V522S29E-05 

18 2.39597te95l3930E-03 2.53O9O20298O-87E-O5 

19 -1.8696569169S250E-03-1. 9605681 76b9<.73E -05 

20 1.98370029882283E-03 1 .59738393B29b97E-05 

21 -l.l85662e9j77639£703-1.231356B2138792E-05 

22 9.985358930B1592E-09 9.61782585020513E-06 

23 -7. 5603195791 BV59E-09-7.80378S-!‘596353E-06 

29 5.97929599892687E-09 6. lS2539250:8793E-06 

25 -9.652565670 120876-09-9.78206561 26821 9E-06 

26 3.590931 15V56986E-09 3.63393633599^05E-0b 

27 -2. 5988992 1222992E-09-2. 6629539 1095'" 39E-06 

28 1.7959937&008782E-09 1 .83799261 759<.S0E-06 

29 -1 .10661968356252E-09-1 . 1 3l226253279«8E-06 

3« 5.13209970392602£-05 5.291 1 053828 79h2E-0V 

THE MATRIX OF FOOMIER COEF. FOR Tnt 6AFFLE COMP. ? 

1 -.66993000763912 - 7.623918621S93V9E-09 

2 -.26663215109372 2.3315209992«.v92E-03 

3 8.82S657e3270970£-02 3.008889897230005-09 
9 -9.6861 26393 1557 lE -02-2. 99682 i29l5Ob90F-0 9 

5 2.962089129dai62E-02 2 . 309a03b999b768E-09 

6 -2.09397038055803E-02-1. 78298e^9023313E-09 

7 1 .56393908729662E-02 1 .91095929281979E-09 

8 -1. 2201 019 1993286£-02' 1.199 1 790bl62771E-09 

9 9.B32l520697l9lt£-03 9.9889621 7262798E-05 

10 -8.130l2S182v2072£-03-8.02320923329V67f-C5 

11 6.871S3021l'bb929£-03 6.90t6260d967379£-05 

12 -5.926315l310a559£-03-6.05151277189395E-05 

13 5-22PC6099389227E-03 5.907b8933512022E-05 
19 -9.72228/95569365E-03-9.96181907291181E-05 

15 4.97718760309763E-03 9. 77978219921 • 19E-0S 

16 -5.07390708385765t-03-5.562l02752S1693E-05 

17 3.19512fc7920l223£-03 3.91996929529552E-0S 

18 -2. 39S97629523659E -03-2. 530902829785595-05 

19 1.86965691836170E-03 I .9605881 766b21 3E-C5 

20 -1 .963?0029890589t-o3-l .59738393b22858£-05 

21 1.18566289669997E-03 I .23135682l90125F-05 

22 -9. 9a535e93i<.639B£-09-9. 81 7825850071615-06 

23 7.5603l9573b2975£-09 7.8o378596557890E-06 
29 -5.97929599993938£-09-6. 15253929998739E-06 

25 9.652565869b6627£-09 9. 7820656127e758E-0b 

26 -3.590931 15993570€-0--3.633936J35922?7E-06 

27 2.5988992ll93b99t-09 2.6b29539l051b07E-06 

28 -I .79599376031 J36£ -09-1 .83799261 750 339r -06 

29 J.1066l96839j592t-D9 I . 1 3122625330V07E -06 

30 -5. 1 32099709 1 9586£ -05-5. 29 11 0538? 73J96E -07 

COMR. input .911352 nozzle EX1.= 5.733113E-03 OISSI.= 
DECAY RAfE* 3.9309577i01-77e5E-C? 


.780698 


decay In decibels /Cycles 2.281 JS72875101 



pROGFiA^' baffle ( INPUT, OUTPUT. TAPE 5 =IKPUTfTAPEf.- 0 UTPUT» 

COMPLEX OMEG,AMU,8MU,ANN.ALAMI,CI,0,ALAMN,CBS,CB1 ,CB?,ECB1 ,ECB2. 
lCCC,esV,BASEv,CCS»BBl.BBt',E81.EB2,BCC,BASEl,Gl .OMEGl » ALIMJ.ALMOZ. 
ZBASE,SLMX,SUMG,SUMF.eSK»A81.AB2.CM.8K 
complex bases 
OIMEKSION TT I (10,10,12) 

DIMENSION 0(A) 

DIMENSION TTE(10,10il2) ,BN ( 1C, ) 0) , AKC ( 1 0, 1 0) ,RN1 U 0 , I 0,90) 
DIMENSION m(M),Z(11) 

OIHENSIPN wl (31), ZI (31) 

COMM 0 N/BLKA/MC,LC,MB,Le,ZB,ALENGTP,AMU(lC, 10 , 12 ) ,BMU( 10 , 10 ) 
COMMCN/BLNO/OISr V,OISF T.T, GAMMA, TM AX 


COMMCN/eLKf/«BES(10,iO),CBES(10,10) 

DATA(tf m , . = 1, 1 n /, 02 129101 e37b5A,.oS76166S831 12 a,. 09 3A0039827825, 
1. 120S2016961 A32, . 1 36A2A90095628, . 1 A 1A9370S92875, . 1 36A2A90095628. 
2.12052016961A32,.C93a 0039827825,.0576I665831 12A, 0212910 18375SA/ 
OATA(Z(I) ,I = l,ll)/.0079S73199523e,.OA69100.77030t>7,.12291663671AS8, 
1.239765 3AA9A7 16,. 360 13A793A 191 1,.5000 o000000000..639815206S803*«, 
^.76923A56550528A,.37/083363285A2,.95308992296933,.9920A2600A7a 2/ 
OATA(wl (I) ,I=l,16)/0.0026887393,0.0075039737,0.012730A?37, / 

10.0176?3180A,0.02?29ao 7S7,0.0267AC7623,0.03100A7839,0.03A92|-96» 
20.038A2AfiAOA,0.0*.15A025lA,O.OAA2e2221S,0.9A6S632991 ,0.?A8321 >'i\ 

' 30.0A95e6799A,0.05038A9228,O.OS06650035/ 

DATA(Z1 (I) ,I=l,l6*/0.00099«8SC7,0.006n037AlO,O.Ol6l30A622, 
10.0313633036,O.OS13677338,0.0758967083*O.IOA790 7a93,0.13?791j 3' 
29.17A5C1629A,0.21a5139136,0.257A590682*0.30292A3265,0.350a0999 Sa, 
30.399AC29S30 ,O.Aa9A289665,0.5/ 

00 78 1=1,15 
Wl(I»16)=V«l (16-1) 

78 Zl (1*10=1. -Zl (16-1) 

the variable MX OETtHM'NcS IE A STANDING WAVE OR TRAVE' ING WAVE 
€•••• SOLUTICN IS 10 BE USEJ. MX=0 WILL GIVE THE STANDING WA '£ AND MX= 
(;*«•• I WILL GIVE THE TRAVELING WAVE **«»****o*ofl****««»****o*«***»**< *« 
MJ;=1 
HX=0 

C««»* INPUTS INCLUDE NO. OF FO'jR IER-BESSEL SERIES TERMS-MC,LC,MG ANO LB. 
L»**« COK2USTOR L£NGT‘' TO RAOI'JS-ALENGTH .EArf LE BLADE LE:-GT:i-ze .BAFFLE 
C**»« blade thicnness-t.nu. or baffle comfartments-mub, ratio of specific 
C o*«« HEATS'CAMHA.WAVE AMPLiTUCE-EPSILON.yEAN FLOW MACS NO.-AMACH, THANS- 
c***« VERSE WOOF CHARACTQR-MHAT. RADIAL MODE CHARACTOR-LHAT . AND MAMMUM 
C*«*« NO. OF ITERATIONS OF THE oUCCESSiVE APPROX. OF SERIES -IDMAX 
READ (5 ,200 » MC,LC,MB,Le, I0MAX,HHAT ,LHAT 

200 FORMAT (7I5J 

READ(5,201 ! ALENGTH,AMACH, GAMMA, EPSILON 

201 FuRMAT (4F 10.5J 
RLA0(5*202» MU8,28,T 

202 FCRMAT(I5,2F 10.5) 

WRITE (6,120) 

120 rOHMAT(*lThE FOLLOWING CALCULATIONS ARE PURELY FOR DIAGNOSTIC 

IFURPOSES*./.* LARGER MATRIX SIZE SHOULu PE USED FOR COMr L‘ AT loN* ) 
WR1TE(6,A0S> 

WKI TE(6,A05) 


I>(MX.£O.OJ A07.A08 
a07 write <6. 909) 

GOTO «:0 

AOe wR!K:6.9n> 

930 C' .' 'InlE 

K . i TE (6 1 9 05 ) 

WRI TE (6 .905) 

909 format (*othe followin:- CALC.sr'f 
911 format (•OThE following calc. Ap. 



FOP STANDING WAVES*) 
FOR traveling waves*) 



AOS FORMAT (»0®.75<1H*) ) 

THESE FaRAMpTERS ARE FOR THE DISSIPATION CALC. 

CTURR=0.05 
VISC=1 .E-05 
OISFv=SC«T»VISC> 

OISFT=SORT (CT'JRB) 

CALL ROOT (MC.LC.2tCBES) 

TALL RC0T(MB,LH.mU3.B3ES» 

c«*** unsteacy combustion Input parameters n and tau are oefineo 

OR=H»-AT 

RLHOA=reE^ (hhAT*j ,L»AT) 

AN= IGAkRA. 1 . ) /A. /GAMMA 
TAU=?. lA •S'i'26S359/RLHUA 
Wfi.‘TF(6,555> 

555 FORMAT («0Th.E CYLINDRICAL CHAM'^EH tiAS THE FOLLOWING CONFIG. •) 

WRITE (6 «SS6» ALCNGTH.AMACH, GAMMA 

556 rORMAT{*0 THE LENGTH TO RAOlUS=».F 10.5.3X.* M£AN FLOW MACH NO.=». 
l-^lP.Si’Tt/.* RATIO OF SPECIFIC HEATS = *«F10.5» 

WRITE i«\s57» LHAT.HHAT.f.PSrLON 

557 FORMAT («0ThE WAVE CHARACTQR iS ••13. • RADIAL MOOE*,3X. 13»» TRANSV 
1 MODE*./.' ANO IS OF WAVE AMPLI rUDE«.G21 • 1A» 

WRITEIc.SGO) MuB.ZB.T 

560 FORMAT (*OThE B&FFLE CONSISTS OF *,I5*« COMPARTMENTS WITH LENGTH*. 
IFIO.S./.' ANO BLADE ThICKNESS‘‘.Fl0.5) 

WRn£(6,SS9» Vise 

559 format :* 0 MOLECULAR V ISCOSlTY-*.G21 .lA) 

WRITEt6,13Al CTURB 

:3A formati’oeffective Turbulent vise. coef. =*, 021 . 14 ; 

PI=3.;a159265359 

0R2-MUB 

DO 100 PU=1>MUB 
OR=MU 

Vl=2.*3. 1A159265359«0R/0R2 
V2=2.*3.14159265359*(0R-1.J UR2 
00 100 MSTAH = 1,.M8 
DO 100 M=1.MC 

ei=<MSTAR-l .)*HL8/2.-IM-1.» 

B0=81«ei 

IFieo.LT.l .E-C6) GOTO 101 
82=(MSTAR-1 . > ‘MLe/Z. • (M-’ .) 

8ASI= (S 1N(B1*V1 ) -SIN<B1*V2) )/(2.«Gl > 

BAS2= IS IN (B2*V1 1-S1N(B2*V3) )/(2.’'B2) 

8AS3=(COS(B. »V1 )-COS(81*V2) )/(2.«Bl) 

BASA ;CCS(B2*V1 )-C0S(B2*V2) )/(2.«C2) 

1F(M>..;0.0) HAs 3 = 3ASA = 0.0 
^TE (KSTAR.M.MU) =BAS1 *BAS2 
TTI (MSTAP.M.MUI =SAS3-RASA 
GOTO UO 
lOi EP=I. 

IF(M.EC.I) EP=2. 

TTE tMSTAR.H.MU) =EP*3. 1A159265359/MU8 
TTI iMST AR. m.mu) =0.0 
100 CONTINLE 
EP = !. 

00 192 ,MB 

0R=IM-1 . ) ®Mi|0/2. 

DO 103 L = ! .LB 
81=381 S (H.Li -• I .E-IO 
B2=BESiCAL (OR.fll » 

8N(M,L>=PI«EP»B2«B2*(1 .-0R*0R/B1/S1 )/MUB 
103 CONTINLE 



10< EP=0.5 
EP=1. 

00 843 Msl«MC 
0R=H-1. 

00 «SS L=l.LC 
e»=CBE£(M,L>*I,E-10 

82=8ESSCAL<0R.«1» 

AHCtPtLl = Pl * EP » e 2 * 82 * Cl «- 0 R * 0 R / Bl/Bll 

155 rChTiHEE 
IFimx.EC.QI EP=0.5 

843 CONTINCE 

THE FCLLOkINC- MATRIX OE INTEGRALS OF PRODUCTS OF BESSEL FUNCTIONS 
C<**« represent a major amount of COMPUTING TIME- FOR A PARTICULAR BAFFLE 

configuration these values neeo be co^puteo only once and storeo 

NQ=MC«LC 
00 1S6 N=1.H0 
00 156 M=1.MB 
00 156 L=1.LB 

156 BN1{M,l.N)=0.0 
00 70 K=1.MC 
MG=H-1 
0RC=M6 

00 72 HP=1,H0 
0R8=tMF-l .)«Mu8/2. 

00 72 L=1«LC 
RLM1=CEES(M»L) 

AL5=3ESSCAL(0RC,PLM1/2.I 
NC=L»LC*hG 
00 72 LP=1 .LB 
RLH2=BEES(MF.LP) 

AL6=8ESSCAL (0K8,RLM2/2. ) 

SUM=0.0 
00 74 1=1.31 
ALX=71 (I) 

ALX 1 = ALX '» PLMI » 2 . 

ALX 2 = ALx ^ PLf '?* 2 . 

IFlORC.GE.S.A'.O.ALXl.LT.OPCl GOTO 300 
AL3=BESSCAL{CRC.RLM1»ALX» 

GOTO 3C1 

300 AL3=AL5* I <2*Alx J «*ORC) 

301 CCNTINLE 

IFIORB.GE. 5. AND. ALX2.lt. ORB) GOTO 302 


CI=ChPLx (O.C, 1 .0) 

00 640 M0=I.MUB 
00 840 H=1 .mb 
00 840 L= J .lb 

840 4KJ(M,L,MU)=(TTE IH,MS.MU)»CI*TTI(M,MS.mu) )*i3N1 (M.u.NC)/ 


AL4 = BESSt AL (0 kB.RLK2*ALX) 
GOTO 303 

302 AL4-- *L6«( c2.«alx)*«0HB) 

303 CONTINUE 

74 SUM=SUH.ALX*m1 (I)»AL3*AL4 
BNI (MP.LP.NC) =SUM 
n CONTINUE 
10 CONTINUE 
10=0 

HS=M»-AT .1 
NC=LHAT.LC*M«AT 
00 845 M=i ,MC 
00 845 L=1 .LC 

845 RMu(M.L}=CPPUX(O.OtO.O) 





c»**« FIRST CUESS OF FREQUENCY 
OHEG=Cr'PLX(l. 0.0.0) 

831 CONTINUE 

NC=LI-AT«LC*NHAT 

OEU-0.1 

HS=Ht-AT»l 

?I=3.I<•159^6S359 

OI!l=Ct'PLX(0.0.-I.O) 

Ole) =CFPLX(-1 .0.0.0) 

QI3J ^Ct-PEXIl .0.0.0) 

CI4,)=C)'PLX(0.0.I.0> 

22 ? TN=0 
215 cO:. IINUE 

ANN=£^* (1 .-CExP(-CI»UHEG«T*U) )-l./OAHMA 

ALXM1 = -GA'->!A*CI*0 *I£G*;.*«ACM*ANN/II..CAwmA«AMACH**MACH*ANN) 

ab=amach» cgawma-i.: ^2. 

ALAKN = -C I»0MEG**e/ 1 1 . ‘AMAt -t«A8‘ 

ces=cscPT ioneg«cmeg*apac) »«.va.: » i ahach«amach-i . ) • (Rl«oa«rlmoa- 

lOHtG^OPEG) ) 

CBl = CCfEG'AHACH^CBS)/ 1 1 .-AhACH*AHACH) 

C62= (O-EG^AHACH-ceS)/ 1 1 ,-amaCH*AmacH) 

ECBl =CE XP ccl* ICB1-CB2) * IZB-ALENGTH) » 

ECF2=CPPLX (1 .0.0.0) 

CCC=-(C I*CB1-ALA«N)/(CI»CB2-ALAMNI 

esv-ci* »CBl»EC61.Cfi2*CCC*ECB2)/IEC81»CCC*ECB2) 

BtSEv-CNPcx(O.O.O.O) 

00 800 

C»= (V-1 . ) •N<jfl/2. 

00 800 L=l.Lfl 
RLHDA8=paf SIB.L! 

CCS=CSCftT !CMt*G*CBEG«ANACH*AMACH* |AMACh*AHACH-1 .)• (RLHCAH*RLH0A8- 
lOMEG-eCVEG) ) 

881 = (C‘'F_G*AMACH»CCS) / C 1 .-AMACH»AHACN> 

982= ICFEC-'AMACh-CCS)/ J1 .-AHACH*AMAC).) 

E82=CEXF<C I* (932-861 :*ZB) 

FB1=C“=LX (1 .0.0.0) 

BCC=-(C I*Rei-A!.A«.I)/(CI*882-ALAHI» 

9ASE1=0*PLX (O.O.O.O) 

00 eei Pb=i.Hu8 

881 8ASE1=EASE1*AV!U(«.L.HU)*<TT£|M.MS .NU) -Cl*TTI|M,MS .Mg)) 
BAS£l=eASEl*3Nl (K.L.NC? 

SUHX= (eei*E9l .se2«£82«8CC)/(£Bl •8CC*ER2) 
eftS£V=EASEV»BAS£l*CI* (BBl*£8l»6B2«£e2*BCC)/(£Pl »6CC*E82) 

BOO CONTINUE 

6AS£v=EaS£V/AMC(MS.U)-AT) 

Gl=aSV-6AitV 
6=CA6S{G1 ) 

IF(In.!E.O) goto 22<* 

223 Oh£G1=CmEG 
REF=G 
IN = 1 

GOTO 2c5 

2Z<i IfiG.Ui.REF) GOTO 22.1 
iF(IN.UT.A) GOTO 228 
OMEGxC^EGl 
0£U=CEL/2. 

:F(0EU.UT.1.£-0A) goto 310 

IN=1 

GOTO 2J5 
228 1N=1N»1 

22; OHEO = OKFGi ‘OT IN)<*0EL 



GOTO 213 

£•••• JHIS IIEHAIEO VALUE OF FHOUENCY ONLY 'NCLUOES THE HAVE ALTERATION 
C*«*« EFFECTS OF ThE 0AFFLE AND CONSEQUENTLY PREDICTS A GROhTh RATE 
C»«*A energy OISSIPATICN hASNT BEEN INTRODUCED AT THIS POINT IN THE PROG. 
310 vHITE lE.i.lO> 

Alc roRHAr(*o;^ESAiicN*.n.* complex freo.=*. 2 G 21 .iai 
C***« iteration of main CSiANatH COEFS. 

A! INjsALAHI 
ALN0Z=ALAMN 
DO 809 MS1=1.MC 
MSO=VST-l 
0R1=MST-1. 

DO 811 LST=1.LC 

IF (MST.EO.hS.ANC.LST.EO.LHATT GOTO 811 
NCl=LST*LC*tHST-l.) 

PLKCa1=CBES««ST.LST1 

BAS£=OREG«OmEG»amacH«Akach. (AMAC h»AMACH-l ,!• ( RLMDA1 «Rlh0AI- 
10MEG«CT'EG» 

BASE=C£CRT <0ASEI 

C81 = C0?'FG»-MACh.BASE)/(1.-AMACH*AMACM1 
C32- :0KEG»AMACh-9ASE)/(1 .-AMACH*AMACH) 

CCC=- CC I«C81-6LN0ZI/ lCI»CR2-ALNO?l 
ECBl =CE aP ( C I • ( C6 1 -C02 J • ( ZB-ALENGTB J » 

ECB2=CPPLX{1.0.0.0J 

CBS=CI« (C61*FCei.C82*CCC»ECB2)/{EC01*rCC*EC82) 

8ASEG=AMCiMST.LST) 

SUMA^CPPLA (O-O.O.O) 

DO 81A M=!.M8 
0R=(M-1.>«muP/2. 

DO 81<i L = 1 .'-0 
RLH=PB£S(M.L) 

CCS=CSCRT (CMEG*CmEG*A«ACm*AMACH. ( AHACh*AHACH- 1 . ) • (RLM*RLM- 
lOKEC'CVEG? > 

8Bl = (0REG*'AMACH*CCS>/n .-tMACH*AMACH) 

8B2= (CyEG^AMSCh-CCSI/ tl .- amach*ahacmi 
EB1=CmFlx(1 .0.0.0) 

EB*'=CEXP(CI* IP02-«B1)*Z8) 

BCC=- (C l»Ffll-ALlNJ)Z(CI*BP2-ALlNJ) 

BASE=C1« .082*FB2*‘»CO/iEB1 »9CC«E82> 

SUHG = C»'FLX(0. 0.0.0) 

DO 817 MU=1.MUB 

817 SU«G=SLRG.aHUiM,L.MUl«(rTE‘M,MST.HU)-Cl«TTl (M,MG. ,HO) ) 

SURG = SLVG»«N1 (M.L.NCH 
819 SUMx=SLMX.SUMO«eASE 

BMU(KST .LST ) =f.uHx/Ui l'/CPS 
811 CONIInLE 
80S COfiTINLE 

PHU(PS.LHAT) =C4FLXU .O.O.f 5 
IO=IC»l 

00 960 H=i.H9 
OR=(K -l.)«PbP/2. 

00 960 L=l.LP 
00 970 PU=1 .PUR 
SUMF=C7’PLX (0.0 •C.')) 

00 971 M.p=l.v: 

00 971 LP=!.LC 
NC=LF*LC« (MP-1 . ) 

SUHF-SLMF •HMU(MP.LP)*BNU«»L.NC)*(TTE(M.BP ,MU)»CI*TTI (H.MP .HU)) 
971 COMINlE 

'WO -SUMF/PN (H.L) 

960 CONTINLE 


^,v; IV 

oV' 


■,0B 



IF(IC.LT.IOMAx) goto 8S1 
WRITE 

909 FORMAT («0Th£ FOURIER-BESSEL COEFS. FOR MAIN CHAMBER** 

DO 911 M=1,HC 

WRITE (( *910) H, (BMUiM.U 

911 CONTINUE 

910 FORMAT!* •.I3,* • * lOGi 0. !./• 10610. 3) 

DO 912 10=1. •‘L'B 

wRirEte. 9*3)10 

913 FORMAT (*0THE FOORIER-8ESSEL COEFS. FOP BAFFLE COMPARTMENT*. I3* 
00 91A M=l,MB 

WRITE (£.910) M« (AHU(M.LtIO).L=l*LB) 

91A CONTINLE 

912 CONTINUE 

c**«* stability calc. 

0MC=REAL(0MEG* 

BlNl=C.0 
SUHV1=0.0 
DO 821 M=I.MB 
0R=(M-1.»*HU8/2. 

DO 821 L=1.LB 
RLK=BeES(M«U 

BSK=CSCRT (0MEG*0MEG*AMACH*AMACH. I AMACH*AMACH-I . * * (RLM*RLM- 
10M€G*0»'EG) ) 

ABlslOMEG'AMACH.eSfO/n .-AMACH*AMAChl 
A82= (CMEG«AMACH'aSM) / I 1 .-AMACH*AMACH» 
CM=-ICI*A81-AL1?«JI/<C1*A82-ALINJ* 

8AS£=CEAP(CI*A81*2B) .CH^CEXP (CI*A82*ZP) 

BASsCAESlBASE) 

BAS=eA«*BAS 

881=!1 .•CM5*C!»CMEG-AHACH*CI*!A8l.Ae2*CM) 

BASR=CA0S«B91 I 

BASR=eASfi'’BASft*GAMMA*GAHMA 

SUMWl=0.0 

DO 822 1=1.11 

ETA=ze*zm 

8SM=CEAP(CI*A81*FTm) .CM^CEaPJCI^ABZ^ETA* 

BS1=CAES{8SM) 

B82=CI* (Aei*CExP(Cl*A61*CTA».CH*A0Z*CEXPIC!*A82*ETA) 

BSM=GAKma« (CI*0MEG*8SM»AHACH*8e2) 

8S =CAES(8SM 

BS=BS*ES/?./GAHPA»GAMMA«0HG*0MG*B51*BS1/2. 

822 SUMwl=S0MWl.BS**Ul)*Z8/8AS 
RAH=0.0 

DO 823 MU=1.MUB 
RAM1=CA9S(A.MU(M,L.MU)) 

823 RAM=KA)'<PAM1*PaHI 

6lNl=elM ♦«A‘<*PASR*8 N(M,l)/8AS 
821 SUKYI=SGmy1 »RAM*sUMw1*8N(N,L) 

8NZI=C.0 
SU*<U1=0.0 
00 841 M=1,HC 
OR=K-l . 

DO 841 L=1.LC 
RLH=C?ES(M,L) 

8K=CSCRT (0M£G*0P£G 'AMACh»A“ACM. |ANAC't«AMACH-l . ) • (RLM«RLH- 
JOHtG«C)'EG) J 

CB1= (CfEG'AMACH.BK)/ ! 1 »-AMACH*AMACH) 

CB2= (OI-EG* AMACH-eft ) / ! 1 .-AmaCH*AMACH> 

CM = - ICI»CB1-ALAHN) / (OC82-ALAMN) 

HASE=CExP (CI» (CBI-C82) • (ZB-ALENGTH) ) ‘CM 



B*S=C*eSfRASE) 

8AS=EAS»BAS 

8AS£e=CExP (-C I«CB2* 128-ALtNGTHl I 
BBS=CAES(RASE8» 

BBS=8eS»6BS 

881 = « 1 . *CM) *C I*CMEG» AMACH*C1* (CBl •C8*CS2» 

8ASR=CABS(881) 

basr=easr*basr«gamha*gahha 

SUHVl=C.O 

00 S <«2 1 = 1.11 

ETA=ze» IALENGTH-2B)«Zm 

0SK=CE>P(CI*(CB1-CB2)*<ETA-ALENGTHI ) *CH 

8S1=CAES(RSA> 

B02=C !• <CRl*CExP (Cl* (CB1-CB2) • (ETA-ALENGTH J I *CM*CB2» 

BSK=CA»'hA* (CI«0PEG*8SK*AMACH*BB2I 
8S =CAES(BSA) 

8S=BS*ES/2./GAMMA*GAMMA»0MG*0MG*BS 1*851/2. 
eASEe=CExP(-CI*C02*(Z8-ETA» ) 
eCS=CAES(BASEBl 
eCS=BCE*BCS 

8A2 SUHVl =S(JMV1 »aCS*PS*H ( I » * ( ALENGTh-28) /8AS 

RAMl=CABS(BBU(M.Ln 

BNZI=6NZI.RAH1*RAMI*BASR*AMC(H,L»*8BS/8AS 
8A1 SU}*U1=SUMU1*RAHI*RAM1oamC(M,L)*SUMV1 
T=T/2. 

CALL VC ISP (0><G. RUB. EPSILON. AMACB. SUM J 

eiMR=ElM/2. 

BNZlfi=eNZI/2. 

FN=A\> ( I .-COS (OMG*TAu) > *AMACH 
ANZB=AKACH*(1 ..GAMMA)/?. /gamma 
01N1R=FN*PIMR 
‘^•.2IR = #NZR*b\ZlR 
SUMIK=EN2 IR-RInIR.SUM 
mR1TE(£.30) BINIR.BNZIR.SUH 

30 FORMAT {*0COMB. INPUT*,G13.6.* NOZZLE EXT.=*.G13.6»* OISSI .=*.G1 3.6 

1> 

SUMY1=SLHU1«SUMY1 

c**** decay rate calc. 

0ECYKT=SUMIN/SUMy 1 
WRITE (t.l?5> OECyRT 
12S rORHAT (*OOECAY RATE=*.G21 .lA) 

OECTe=5A.S75*DECYRT/OMEG 
wRITE(e.l30> 0ECT9 

130 FORMAT (*ODECAr IN DEC IBELS/CYCLE=* «G2 1 . 1 A 
END 



SUBROUT I NE VO I SP ( OMC G . HUB * EPS I LON t A <4 AcH • V ALUE ) 

COMPLEX A,4HUt0HU»SUM,SUMX 
DIMENSION 4 «l?t’ 1 I 
DIMENSION N(ll) <111 

COHHCN/8LKA/MC.LC<HB«LB«ZB«ALCNGTH*AHu<10tlO*l2) tBMU<10*10> 
C0MH0 N'/PLi' 0/0ISEV»0ISFT.T,GAMMA,TMAX 
COMHCN/PLKr/ePES(10*10) «C0ES(IO*1O> 

OATA(h(I) «I = 1 ,1 n/.0ZlZ9101837SS<> *.05761 665831 124. .0934400398278?5t 

1.12052016961432.. 1364?490095628».!414937089?875».136424900956?8. 

2. 12052016961432. . 09340039827825..05761665831 124..021 2P1 01837554/ 
OATAlzm « I = !*I D/.00795731995258..04691007703067..122916636714S8. 

1 .23076534494716. . 360 184 7934 1911 ».5000p000000000» .63981520658089. 

2.769234565505284. . 87708336328542.. 95308992296933. .9920426004742/ 
C««o* DISSIPATION C4LCULATI0N IS AS FOLLOWS 

81=3.14159265359 

c««** smallest value for RA0MIN=Mu8*T/PI 

RADHIN=MUB*T/PI 

DO 331 1=1.11 

R=RACMlN» (1.-RA0 mIN)*Z(H 

TMAX=T/R 

DO 331 MU=1.MUB 

THETAX=2.®PI*MU/FL0AT (MUSI »TMAX 
THETA =2.tpl*Mu/FLOAT(HUB)-TMAX 
Th=2.*Fl«HU/FLOAT (MU8) 

AL=He*L8 

SUMX=CPPLX (0.0.0. 0> 

SUM=CMPLX (O.O.O.OI 
DO 332 H=l.Mfl 
OR= .8-1 . ) «Mue/2. 

DO J32 L=1.L8 
V=BBES(M.L) 

V2=BESSCAL(0R,V/2.> 

V=V*R 

V3=V*2. 

IF(OR.CE.5.ANO.V3.LT.OR) GOTO 300 
V1=BESSCaL (OR.v) 

GOTO 301 

300 Vl=V2*((R*2.>**CR) 

301 CONTINLE 
M0X=MU«1 

IF (MU.EO.HuBl MUX=l 

SUMX=SLHx*AMU (H.L .MUX)»V l»COS(OR*THETAX) 

332 SUH=SUF»AHU(M.L,HU) • VI ‘COS (OR*Th£TA> 

A (HU. 11= (SUMX-SUHl/SQRT (2.*TI 
331 CONTINLE 

AMULT=GAHhA»SQRT (0HEG/2.)/2. 

SUM=CmFLX(J. 0,0.0* 

DO 339 J:-l,ll 
TE=ACOS i2B/(2B.T) » 

DO 339 HU=l.Mufl 
AX=CA8S (A (MU, J1 ) 

AX=Ax*Ax 

B=EPSILCN*EPSIL0N»AX/T/4. 

SUMI=0.0 
DO 336 1=1.11 
ETA=TE/2.*Z ( I)*(Pl-TE)/2, 

V1=SIMETA) 

336 SUMI=SLMI.(DISFV-DISFT*( ( 4MACH»AMACH»P*V1*V1 ' •*0.25» )»W(l)/Vl 
SUMI=SLMI«4X*(PI-TE»/2. 

339 SUH=SUH .SUMl«w ( J) • ( 1 .-RAOMINI 
VALUE=AHUL I*SUM 
RETURN 
END 



FUNCTICN BFSSCAL (ORtALxi 

C**»* this StBPROGRAM CALCULATES THE VALUE OF THE BESSEL FUNCTION Of 
integer and half INTEGER OROER-OR WITH ARGUMENT -ALX 
M=OR 

A3=(OR-m«(OR-M» 

1F(A3»C-T.0«n01) GOTO 100 
IF(ALX.GE.3.0) GOTO 101 
AL2=ALX«ALX/3./3. 

ALA=ALi»AL2 

AL6=AL^*ALA 

AL8=ALA«AL4 

AL10=AL8*AL2 

AL12=AL6»AL6 

BASEl = 1.0-2.24q9<J97*AL2»1.2656208«AL4-0.31638e>6«AL6»0.0<iA<i479«AL8- 
10.003 S<i44.*AL10‘0.0002100«AL12 

8ASE2=O.S-0-5629<}985»AL2.0.21093573»ALA-0.0395A289*AL6»0.004«3315* 
1AL8-0.C0031 761«AL10»0.00001 109»AL12 
BASE2=ALX*BASE2 
GOTO 200 
101 AL1=3./ALX 
AL2=AL1*AL1 
AL3=AL2<AL1 
ALA=AL2*AL2 
AL5=AL2*AL2 
AL6=AL2*AL3 

F 1 =0. 79 788.956-0 > 00000077* AL1-0.005527a*AL2-0. 00009512* AL3* 

10. 001 37237* Al9-0.00072B05*AL5» 0.000 199 76* AL6 
THETA1=ALX-0. 78539816-0. 09166397*AL 1-0. 0000395A*AL2*0. 00262573* 

1 AL3-C.C 00591 25* AL9-0. 00029333* AL5»0. 000 1 3558* AL6 
F2=0.7? 78P9S6*0.000001S6*AL1 *0.016596^7*Al 2*0.0001 7I05*AL3* 
10.0029951 1«AL9. 0.00 11 3653* ALS-0. 00020033* AL6 
THETA2=ALX-2.356l9999*0.l299O6i^'*Ll»0.00305650*AL2-0.00637879*AL3 
1* 0.00079 398* AL9*0. 000 79829*AL5-0. 000291 66* AL6 
B1=ALX**0.5 

BASE1=F4*COS(ThETA1)/81 

BASE2=F2*COSjTHETA2)/81 

200 IF(M.EC.O) GOTO 201 
IF(M.FC.l) 202.203 

201 BESSCAL=BASE1 
RETURN 

202 BESSCAL=BASE2 
RETURN 

203 CONTINLE 
Ml=M-l 

DO 209 K2=l .Ml 
BASE3=2 .*h2*BASE2/ALX-BASE1 
BA5E1=FASE2 
BASE2=E#SE3 
209 CONTINLE 

BESSCAL :3AsE2 

return 

100 CONTINLE 

BASE=(2.0/J.1915926/ALX)**0. 

BASEUCOScALXf^BASE 

BASE2=SIN(ALX)*BAS£ 

IFIH.EC.OI 220.223 
220 BESSCALr8ASE2 
RETURN 
223 H1=M 

00 229 M2=l .Ml 
AM3=M2-0.5 



Bi»SE3=c.*AM3»8ASE2/ALJl-0ASEl 

8ASE1=EASE2 

BASE3=eASE3 

CONiINLE 

BESSCAL=8ASE2 

RETURN 

END 


SUBROUTINE ROOT jMX»LX»MUB.R 
DIMENSION RT(10*10> 

DO 120 M=l,Mx 
0R=MU8«{M-1.)/2. 

OR3=OR**0. 33333 
RT(H,1)=0.0 
lElM.EC.l) GOTO lAO 

RT(M, l)=OR« 0.0086 165*OR3*0.072AS/OR3-O.OS097/OR*0.009A*OR3/''° 
140 CONTINLE 

DO 121 L=2.Lx 

Z=RT (H,l-D ♦3.1415926535 

0IR=1 

OEL=0.1 

REF=CR«BESSCAl (0R.Z)/Z-8CSSCAL(0R*1»Z) 

100 Z=Z»CEL*OIR 

REFl=CR*BESSCAL(ORfZ)/Z-0ESSCAL(OR*l*Z> 

V=REE1/REE 
IF(V.GT.O) GOTO 100 
1F(OEL.LT.1.oE-04) GOTO 105 
OEL=DEL/2. 

dir=-cir 

REF=REF1 
GOTO 100 
105 RT(M,L)=Z 
121 CONTINLE 
120 CONTINLE 
RETURN 
END 



the following CALCdLAI ;oN5 are purely for diagnostic purposes 
larger matrix size S^'CuLO be used for computation 


THE rOLLOwING CALC. ARE FOR STANDING WAVES 


THE Cylindrical chamber has the following config. 

THE length IC RA0IUS= l.SOOOO MEAN FLOW MACH NO.= .10000 
RATIO OF SPECIFIC H£ATS= I .^0000 

THE WAVE CHAhACTOR is 1 radial MODE I TRANSV MODE 

AND is OF WAVE AMPLITlOE .10000000000000 

THE baffle CONSISTS OF 3 COMPARTMENTS WITH LENf 'H .50000 

AND BLADE thickness .lOCCO 

MOLECULAR VlSCOSlTY= 1 • 00000 00000000 OE-05 


EFFECTIVE TURPULENT ViSC. CCEF.= 5.00000000000000E-03 


ITERATION 0 Complex FRE0.= I .666A06as9R99V 

ITERATION 1 complex FRE0.= 1 .58‘50fe2<,999999 

ITERATION 2 Complex FRE0.= 1 .62A02J<.37<.999 

ITERATION 3 COMPLEX FREO.= 1 .608E0312<.9999 

ITERATION <1 Complex FRE0.= 1 .615A29G87<.999 

THE FOURIEP-PESSEL CCEFs. FOR MAIN CHAMBER 

1 -A.0f6E-16 1.052E-15 2.C2AE-1A 1.093E-1«. 

2 1.00 0. .133 1.017E-02 

3 -.2A5 -3.'.SaE-02-9.036E-02-2.309F-03 

A -2.23 3E-13-7.3?9E-1<—5.905E-1A-<..09 7E-15 
5 .121 7.6A1E-03 5.256E-02-1 .791E-03 


-Y.61 32812A999996E-02 
-S.2S39062A9R9991E-02 
-3.9f,A8A37A9'39988E-02 
-A.8632ei2<.999962E-02 
-A.2773A374999980E-02 


THE FOURlER-tiESSEL COEFS. FOR gAFFLE COMdaRTHENT 1 

1 .222 2.9bS£-03-9.Sl<.E-02 S.927E-0A 

2 .540 -3.071E-02 .109 8.315F-03 

3 -.200 -1.045E-02-8.382E-C2-4.hS9E-03 

4 .209 1.123E-02 9.777E-02 6.530E-04 

5 -8. 90 1 E-02 -4. 356E- 03-4. 6 16£- 02- 1 .0S6E -03 

THE FOURIER-c'EssEL CCEFS. FOR BAFFLE COMPARTMENT 2 

1 -.4<,3 -5.9nE-03 .190 -1 .08SE-03 

2 -1.814E-12 4.e3flE-15-3.64lE-14-l .521F-1S 

3 .400 2.091F-02 .lf-8 9.319E-03 

4 -2.7ME-13 8.218E- 15-1 .lSlE-13 9.525E-16 

5 .178 8.712E-03 9.232E-02 2.172E-03 


THE FOURIER-BEsSEL CCEFs. FOP BaFFuE COMPARTMENT 


1 .222 

2 -.540 

3 -.200 

4 -.209 


2.95SE-03-9.S1 aE- 02 5.427F-04 
3.C7u:-02 -.109 -8.315F-03 

-1 .045F-02-8.382E-02-<..f,S9r-0 3 
-!.123E-02-P.777f-0?-6.530E-04 


5 -li.90lE-02-4.35fcF-0 3-4.81hE-02-l .086F-0 3 


3 


COMB. Input .108318 nozzle EXi.= 3.177977E-02 OISSI.= 


OECA. RATE= .2090583E299J36 



.418200 


CE^Ay In OEClEELS/CrCLEr 7.CS77917121102 




